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[57] ABSTRACT 

A radio communication system and method for nnnimizing 
co-channel interference comprises a technique for commu- 
nicating between first stations and second stations. The 
technique involves: i) processing, at a first station, first 
signals into a form representative of modulated radio sig- 
nals; ii) combining the representations to form a number of 
transmit signals using a set of combining parameters; iii) 
transmitting the transmit signals using a transmitter and 
antenna for each of the signals; iv) receiving the transmit 
signals at the second stations to produce received signals; v) 
measuring at the second stations an amount received of at 
least one unwanted signal relative to the amount of wanted 
signal received to produce measurement results; vi) trans- 
mitting the measurement results from the second stations to 
the first station; and vii) receiving the measurement results 
at the first station and using the results to modify the 
combining parameters such that unwanted signals received 
at the second stations are reduced and wanted signals 
enhanced. The measuring step includes the step of correlat- 
ing received signals with known signals to produce the 
amount of at least one unwanted signal, and the amount of 
wanted signal. 

22 Claims, 39 Drawing Sheets 
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CELLULAR/SATELLITE transmitters and antennas for three cells are collected onto 

COMMUNICATIONS SYSTEMS WITH toe same site for economy, the antenna patterns are then only 

IMPROVED FREQUENCY RE-USE required to illuminate 120 degree sectors, and the resultant 

azimuthal directive gain largely compensates for the double 

This application is a continuation of application Ser. No. 5 distance to the far side of the cell. The antenna pattern can 

08/442,648, filed May 16, 1995, now abandoned, which, in be shaped appropriately so as to provide a gain commensu- 

turn, is a divisional of application Ser. No. 08/179,953, filed rate with the maximum range needed in each direction, 

Jan. 11, 1994, now U.S. Pat. No. 5,619,503. which is halved at +/-60 degrees compared to mid-sector. 

Thus a sectorized antenna pattern can be narrowed to -12 dB 

BACKGROUND 10 at +/-60 degrees, giving a mid-sector gain of about 8 to 9 dB 

The present invention relates to radio communication to ^ achieving the maximum range in that direction, 

systems with increased capacity. The system can include a Using central illumination, the U.S AMPS cellular mobile 

number of roving, automobile-mounted or handheld tele- telephone system denies re-use of the same frequency within 

phone sets served by either fixed, ground-based stations or a 21-cell area around a given cell. This is called a 21 -cell 

by orbiting satellites or by a combination of both. The 15 frequency re -use pattern and results in co-channel interfer- 

capacity of such systems to serve a large number of sub- ence being approximately 18 dB below a wanted signal 

scribers depends on how much of the radio spectrum is wnea aU channels are concurrently in use (commonly called 

allocated for the service and how efficiently it is used. maximum load). Such a 21-cell re-use pattern is illustrated 

Efficiency of spectral utilization is measured in units of ^ ^G. 2. Certain re-use pattern sizes such as 3, 4, 7 and 

simultaneous conversations (erlangs) per megahertz per 20 products thereof (e.g., 9, 12, 21 ... ) result in co-channel 

square kilometer. In general, spectral efficiency can be interferers being equidistant from the wanted signal and 

improved more by finding ways to re-use the available located on the vertices of a hexagon, separated by a number 

bandwidth many times over than by attempting to pack more of cells equal to the square root of the pattern size, 

conversations into the same bandwidth, since narrowing the In practice, illumination takes place from sites at the 

bandwidth generally results in the need to increase spatial 25 junction of three cells. Although the re-use pattern is a 

separation between conversations thus negating the gain in 21-cell pattern, it can also be described as 7-sites each 

capacity. Therefore, it is generally better to increase the having a 3-frequency re-use pattern around the three, 120 

bandwidth used for each conversation so that closer fre- degree sectors. The signal to co-channel interference char- 

quency re-use is possible. acteristics arising from this form of illumination are not 

Spread-spectrum communications systems (e.g., CDMA 30 exactly equivalent to those characteristics which result from 

systems) that increase the signal bandwidths using heavy central illumination (due to the antenna directivity it can be 

redundant coding, such that a signal can be read even shown that interference with respect to a particular signal 

through interference from other users, offer high spectral arises principally from two other sites whose antennas are 

efficiency. Using such systems, several users in the same cell fi rin g m me right direction, and not from six equidistant cells 

can coexist in the same bandwidth, overlapping in both which transmit on a common frequency as would be the case 

frequency and time. If co-frequency interferers in the same m central illumination). 

cell can be tolerated, co-frequency interferers one or more The 3-sector, 7-site method of illumination is sometimes 

cells away can also be tolerated since distance will lessen called sectorization, which can give the erroneous impres- 

their interference contribution, so it would be possible to sion that an originally larger cell was split into three smaller 

re-use all frequencies in all cells. cells or sectors by use of directional antennas. This 

Spread-spectrum system capacity is said to be self- impression, however, is inaccurate because the arrangement 

interference limited because each unwanted signal that is used for illuminating three cells from the same site is merely 

received simultaneously with the desired signal, and on the an economic arrangement that actually has slight disadvan- 

same frequency, contributes an interference component. 45 tages over central illumination with respect to technical 

Some systems, however, such as satellite communications performance but is otherwise very similar, 

systems, are already limited by natural noise, so the wide- Cell-splitting is another concept entirely, being a way of 

band spread-spectrum approach is then not necessarily the obtaining more capacity per square kilometer by providing 

best technique for maximizing capacity. Consequently it base stations more densely on the ground. Introducing cell 

would be desirable to re-use the whole spectrum in every 50 splitting in an already existing system usually requires 

adjacent cell or region without incurring the self- complete revamping of the frequency re-use plan, as it is 

interference penalty of wideband spread-spectrum. conventionally not possible simply to split a cell, for 

FIG. 1 shows a typical arrangement of a cellular telephone example, into three cells and to re-use the original frequen- 

network using land-based stations. This figure is illustrative cies three times over. This would result in the three new cells 

of such networks only, for example, cells are not always of 55 operating on the same frequency with no spatial separation, 

such regular size and shape and as a general definition a cell which would present a problem for a mobile phone on the 

may be described as an area illuminated with a distinct boundary between two cells where it receives equal strength 

signal. (but different content) signals on the same frequency from 

Cells can be illuminated from their geographical centers, DOta - it would be desirable to allow a cell to be split 

but it is more common to illuminate a cluster of three cells 60 sectors with the same frequencies being used in each 

from a common site at the junction of the three cells, as site without the above-described interference problem, 

real estate cost is a major economic consideration. The Similar capacity issues arise in designing a satellite com- \ 

antenna radiation patterns for central illumination of a cell munications system to serve mobile or handheld phones. On 

would generally be omnidirectional in azimuth. It is also handheld phones, omnidirectional antennas of indifferent 

common to narrow the radiation pattern in the vertical plane 65 performance are all that in practice the majority of consum- 

so as to concentrate the energy towards land -based tele- ers are willing to accept. Directional antennas that have to be 

phones and avoid wasting energy skywards. When the oriented toward the satellite or larger, more cumbersome 
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antennas do not now find favor in the marketplace, so it is 
necessary for the satellite to provide a high enough signal 
strength at the ground to communicate with such devices. 
The signal strength received at the ground from a satellite is 
usually measured in units of watts per square meter or dBW 
per square meter on a logarithmic scale. For example, a flux 
density of the order on -123 dBW per square meter is used 
for voice communication to provide an adequate link margin 
for multipath fading, shadowing, polarization mismatch etc., 
using a downlink frequency of 2 GHz. The total number of 
watts radiated by the satellite is then equal to this required 
flux density times the area of the geographical region it 
illuminates. For example, to provide such a voice channel 
anywhere in the entire United States, having an area of 9 
million square kilometers requires a total radiated power of: 

10 _12 - 3 x9xl0 12 =4.S watts from the satellite. 

One voice channel would not, of course, provide a useful 
capacity. Five to ten thousand Erlangs is a more reasonable 
target for serving the United States. One way of increasing 
the capacity would be to generate 4.5 watts on other fre- 
quencies too, each of which could carry one voice channel; 
but a 45 k watt satellite would be very large and expensive 
to launch and would not be an economic way to provide 
10000 erlangs capacity. It is therefore more efficient, having 
used 4.5 watts of satellite RF power to create one voice 
channel's worth of flux density at all places in the United 
States, to find ways which will allow the voice carried by 
that flux to be different at different places, thus supporting 
many different conversations using no more power or band- 
width. 

The ability of a satellite to modulate the same radiated 
flux density differently in different directions depends on the 
angular discrimination provided by its antenna aperture. The 
angular discrimination of an antenna (in radians) is on the 
order of the ratio of the wavelength to the diameter of the 
antenna. Using an exemplary downlink frequency of 2 GHz 
(15 cm wavelength) an antenna of 1.5 meters in diameter 
theoretically has an angular discrimination on the order of 
Vioth of a radian or 5.7 degrees, which, from an orbital 
height of, for example, 10000 kilometers, allows discrimi- 
nation between 37 different directions within the United 
States coverage area. Thus, the same 45 watts of satellite 
radiated power could then support not just one, but 37 
different conversations. 

One way of creating 37 different beams is shown in FIG. 
3. A parabolic reflector focuses the radio energy from a 
pattern of 37 different feeds down to the earth. An image of 
the feeds is projected onto the ground forming the desired 
separately illuminated areas. Unfortunately, using this tech- 
nique there is spillover from one area to another, and in any 
case a mobile phone on the boundary between two or three 
cells receives equal signals from two or three feeds. If these 
signals are independently modulated, the phone receives a 
jumble of three conversations which it cannot decipher. 
Accordingly, conventional systems have been unable to 
exploit the potential capacity increases which would be 
realized using discrimination. 

SUMMARY 

These and other drawbacks and difficulties found in 
conventional radio communication systems, satellite com- 
munication systems and hybrids thereof are overcome 
according to the present invention. 

According to exemplary embodiments of the present 
invention, matrix processing can be used to form numerical 
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combinations of data sample streams. The matrix coeffi- 
cients are selected, and can be periodically adjusted, so that 
each of a plurality of receivers receives its intended signal 
with substantially zero interference. 

5 According to another exemplary embodiment of the 
present invention, signal processing does not adapt to the 
movement of mobile phones or to new call set-up and 
termination, but operates in a deterministic way and instead 
the traffic is adapted to the deterministic characteristics of 

10 the signal processing using a dynamic traffic channel assign- 
ment algorithm. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing, and other, objects, features and advantages 
15 of the present invention will be more readily understood 
upon reading the following detailed description in conjunc- 
tion with the drawings in which: 

FIG. 1 illustrates a conventional land based cellular 
network; 

20 

FIG. 2 illustrates a conventional 21-cell frequency re -use 
plan; 

FIG. 3 shows a conventional satellite implementation of 
37 beams illuminating a region of the Earth; 
25 FIG. 4 illustrates an illumination pattern used to describe 
a feature of the present invention; 

FIG. 5 shows a 3 -cell frequency reuse plan; 
FIG. 6 shows a satellite -mobile communication system 
according to an exemplary embodiment of the present 
30 invention; 

FIG. 7 illustrates a mobile-to-hub transponder according 
to an exemplary embodiment of the present invention; 

FIG. 8(a) illustrates a hub-to-mobile satellite transponder 
3S according to an exemplary embodiment of the present 
invention; 

FIG. 8(6) illustrates a combining network for a power 
amplifier matrix according to another exemplary embodi- 
ment of the present invention; 
40 FIG. 9 shows a hubstation according to an exemplary 
FDMA embodiment of the present invention; 

FIG. 10 illustrates coherent beam signal transmission 
according to an exemplary embodiment of the present 
invention; 

45 FIG. 11 shows spectral characteristics using dual polar- 
izations on k-bank hublinks according to an exemplary 
embodiment; 

FIG. 12 is a block diagram illustrating phase coherent 
transportation of the beam signals according to an exem- 
50 plary embodiment; 

FIG. 13 is a block diagram illustrating phase coherent 
transportation of beam signals according to another exem- 
plary embodiment of the present invention; 
5S FIG. 14 illustrates mapping of 2-bit multiplexed I and Q 
signals to a K-band carrier vector; 

FIG. 15 is a block diagram illustrating yet another exem- 
plary embodiment of phase coherent beam signal transpor- 
tation; 

50 FIG. 16 is a block diagram illustrating hubstation transmit 
signal processing according to an exemplary TDMA 
embodiment of the present invention; 

FIG. 17 illustrates connections between a receive control 
processor and a transmit control processor according to an 
65 exemplary embodiment of the present invention; 

FIG. 18 shows a land-cellular exemplary embodiment of 
the present invention; 
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FIG. 19 is a block diagram illustrating the maximum 3. FIG. 4 illustrates a cross section of the illumination 

likelihood demodulation of signals from an antenna array intensity produced on the ground from an antenna, such as 

according to an exemplary embodiment of the present the antenna shown in FIG. 3. Even for a best interference 

invention; case where a mobile unit is located at the center of beam 2 

FIG. 20 shows an exemplary arrangement of staggered 5 (point A) the illumination from beams 1 and 3 is not zero, 

sector patterns; Dut on ty somewhat reduced. The total signal received by 

FIGS. 21(a) and 21(b) illustrate progressive illumination mo J ilc 2 can * described as the sum of three components, 

patterns according to an exemplary embodiment of the s as * 

present invention; an amount C21 times the beam 1 signal SI (small) 

FIG. 22 is a block diagram illustrating part of an exem- 10 an amount C22 times the beam 2 signal S2 (large) 

plary implementation of a dynamic channel assignment an amount C23 times the beam 3 signal S3 (small) 

embodiment of the present invention; Considering now the reverse (uplink) direction, and 

FIG. 23 is a graphical representation of an exemplary assuming reciprocal propagation, the satellite receives in 

radiation pattern for a circular-symmetric, uniform aperture 1S Deam 2 a contribution from three mobiles, namely C21.M1+ 

illumination function; C22.M2+C23.M3 where Ml, M2, M3 are respectively the 

FIG. 24 is an exe'mplaiy graph of relative signal gain ^? nals radiated from mobiles in cells 1, 2 and 3. If mobile 

versus beam crossover points; 1 15 not close t0 the of beam 2 ' 021 wdl be sma ^ smce 

_ _ „ . , , ... . _„ mobile 2 is within beam 2. C22 will be large; and if mobile 

FIG. 25 is an exemplary graph illustrating C/I versus 3 is not close to the edge of beam 2, C23wiU be smaU. Thus 

mobile position m cell for a 3-cell frequency re-use pattern; 20 ^ ^ ^ ^ ^ ^ ^ ^ ^ 

FIG. 26 is an exemplary graph illustrating C/I versus of cells, it may be that the level of inter-ceU interference can 

beam edge crossover point; be tolerated. 

FIG. 27 is an exemplary graph illustrating C/I versus On the other hand, if a mobile unit is, for example, close 

mobile position in cell for an immediate frequency re-use 2S to the boundary between cell 1 and cell 2, the coefficient C21 

system; will be large and Ml will interfere with the decoding of 

FIG. 28 is an exemplary graph illustrating C/I versus signal M2. The conventional technique for avoiding this 

beam edge crossover point for an immediate frequency problem is to deny use of the same frequency in immediately 

re-use system; adjacent cells. For example, the 3-cell frequency re-use 

FIG. 29 illustrates an exemplary radiation pattern for a 30 £" era shc T in FIG * 5 ^ ^ shaded ™ 

circular-symmetric, * cosine aperture illumination function; ™- 5 ™ ^ft* 51 ^ a fi J st fr^ncy H, while the other 

™„ ~ A . , . - , . .,. cells use £2 and £3 in the indicated pattern. It can be seen that 

FIG. 30 is an exemplary graph of relative signal gam ^ ^ ^ same fr fl do not abut and have 

versus beam crossover points for the dlumination function edg^.^ separations of just less than one cell diameter, 

pattern o , ^ mobile on the edge of one beam is quite far down the 

FIG. 31 is an exemplary graph illustrating C/I versus illumination intensity curve of other cells using the same 
mobile position in cell for a 3-cell re-use pattern for the frequency, and thus avoids interference. However, the draw- 
illumination function of FIG. 29; back is that only one-third of the available frequencies may 

FIG. 32 is an exemplary graph illustrating C/I versus be used in each cell, reducing the spectral utilization effi- 

beam crossover point at all mobile positions within 25% of 4Q ciency by a factor of three. Accordingly, the present inven- 

cell radius for a 3-cell re-use pattern for the illumination tion provides, among other features, a means of cancelling 

function of FIG. 29; the co-channel interference without the loss of spectral 

FIG. 33 is an exemplary graph illustrating C/I versus efficiency entailed in denying spectral re-use in adjacent 

mobile position in cell for immediate frequency re-use for cells. 

the illumination function of FIG. 29; 45 If the expressions for the signals received in all beams Bl , 

FIG. 34 is an exemplary graph illustrating C/I at all points ^2, B3 . . . etc., are collected together and we assume for 

within 25% of beam radius as a function of dB down at beam * c ™ m ?" thc ^ numbcr of mobllc 1 sl S nals ™ ^zms, 

edge crossover for an immediate frequency re-use system ^en the following set of equations results: 
using the aperture illumination function of FIG, 29; 

FIG. 35 illustrates beam and cell patterns according to an 50 
exemplary embodiment of the present invention; 

FIG. 36 illustrates another exemplary beam and cell 

pattern using seven communication channels; " [ 

FIG. 37 is a block diagram of a fixed beam forming 55 ' " 

apparatus according to yet another exemplary embodiment Bn ~ Cnl ' Ml + Cn ^ • A/2 + c«3 • A/3 + cm • ma . . . &m • Mn 
of the present invention; 

FIG. 38 is a diagram of current injection and extraction which can be abbreviated to B=C.M where B and M are 

points illustrating a beam forming apparatus according to an column vectors and C is a square, nxn matrix of coefficients, 

exemplary embodiment of the present invention; and 60 From the signals received by the satellite in each of its 

FIG. 39 illustrates an exemplary TDMA embodiment of beams, it is desirable to determine the signals transmitted by 

the beam forming apparatus of FIG. 37. the mobiles; according to the present invention this can be 



Bl = Cll-Afl+Cl2-A/2 + C13-A/3 + C14-A/4...Cl«-A/« 
B2 = C21 • A/1 + C22 • A/2 + C23 • A/3 + C24 • A/4 . . . C2/i • M/i 
53 - C31-AA +C32- A/2 + C33- A/3 + C34- A/4. .. C3n- A/n 



DETAILED DESCRIPTION 



done by solving the above set of equations to obtain: 



Initially, it is helpful to understand the interference prob- 65 
lems associated with the transmission of signals from con- This solution can be obtained as long as the matrix C is 

ventional antenna arrays, such as the one illustrated in FIG. invertible (i.e., has a non-zero determinant) and results in 
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cancellation of substantially all interference between mobile 37x6 or 222 MHz of bandwidth will be needed on the 

signals and complete separation therebetween. All the ele- satellite-hub link. Since one method of maintaining coherent 

merits of the above equations, that is mobile signals Mi, beam signal transport uses at least twice this bare minimum 

beam signals Bk and matrix elements Cki, are complex bandwidth, and the reverse direction requires the same 

numbers having both a real and an imaginary component so 5 amount, 1 GHz of bandwidth is needed. This suggests that 

as to be able to represent not only signal amplitude differ- a carrier frequency around, for example, 20 GHz is appro- 

ences but also signal phase relationships. According to the priate for the satellite-hub forward and return links, 

present invention the signals received in the different At such a frequency, even relatively small hubstation 

antenna beams are sampled at the same time at a rate dishes will have very narrow beamwidths, so that exclusive 

sufficient to capture all signal components of interest accord- 10 use of this bandwidth by any one system is not necessary, 

ing to Nyquist's criteria. One set of such samples forms the and the entire bandwidth can be re-allocated to other satel- 

column vector B at any instant, and each such vector is lites and ground stations without interference as long as the 

multiplied by the inverse of C, for example, once per sample sightline from a first ground station to a first satellite does 

period to obtain a set of samples M representing interference not intersect with a second satellite. This can be avoided by 

free mobile signals. Successive values of the same element 15 allocating unique "stations" to satellites in geostationary 

of M form the sample stream corresponding to one mobile orbit, or, in the case of lower orbiting satellites that move, 

signal. This stream is fed to a digital signal processor for the probability of intersection is low and can be handled by 

each mobile signal that turns the sample stream into, for having an alternative hub location which is activated when 

example, an analog voice waveform or 64 KB PCM digital such an event threatens. 

voice stream as required by the telephone switching system 20 FIG. 7 shows a block diagram of an exemplary satellite 

to which the system is connected. transponder for relaying mobile -originated signals to the 

According to another aspect of the present invention, the hubstation. The L-band (e.g., 1600 MHz) multi-beam satel- 

matrix C does not have to be inverted every sample period, lite antenna 470 receives signals from a plurality of mobile 

but can be inverted less frequently or only once at the phones distributed between the various beams and amplifies 

beginning of a call. The matrix C and its inverse vary 25 them in respective low-noise amplifiers 480. The composite 

relatively slowly because the rate at which the C coefficients signal from each beam contains, for example, signals from 

change due to the mobile unit shifting position within the 400-500 mobile phones using different frequencies spaced 

beams, or due to the beam illumination intensity distribu- at 12.5 KHz intervals over a total bandwidth of 6 MHz. The 

tions changing due to satellite movement in the non- composite signals of each beam are downconverted in 

geostationary case, is relatively low. In an exemplary satel- 30 respective mixers 440 to obtain baseband signals, for 

hte embodiment of the present invention, typical cell sizes example, spanning the range of 1-7 MHz. This type of 

are in the hundreds of kilometers diameter range and satel- signal will be referred to hereafter as a "video" signal as it 

lites orbiting at medium altitudes take an hour or two to pass is typical of the frequency range spanned by signals from a 

by a typical cell. Thus the need to compute a new matrix TV camera. To downconvert the composite received signal 

inverse due to movements may not arise for the duration of, 35 to the video signal, the downconverters can, for example, be 

for example, a typical 3-minute telephone call. The principal image-rejection type downconverters. The downconversion 

reason that changes in the inverse C-matrix would be process can optionally take place in one or more steps using 

beneficial, however, is that conversations are continually appropriate intermediate frequencies, 

being connected and disconnected. If n=37, for example, and The downconvertors in the system can use the same local 

the average call duration is 3 minutes, then on the average 40 oscillator signal so as to preserve the phase relationships at 

one mobile and its corresponding column of matrix C drops the downconverted frequencies that were received at the 

out and is replaced with another column of coefficients every antennas. The inadvertent introduction of fixed phase mis- 

5 seconds. The process whereby new inverse C matrices are matches and small amplitude differences between channels 

introduced when this occurs will be explained later, suffice is not a problem as this can be calibrated out by analog or 

it to say that this represents a relatively negligible compu- 45 digital processing at the hubstation. 

tational effort compared to the total digital signal processing The baseband signals are used to modulate respective 
involved in demodulating and decoding 37 mobile signals. carriers at the satellite-hub frequency band, e.g., 20 GHz. If 
An exemplary embodiment which applies these principles single-sideband modulation of a 1-7 MHz "video" signal 
will now be described with reference to FIGS. 6-11. were applied to a 20 GHz carrier frequency, the resulting 
a FIG. 6 illustrates a plurality of portable stations 420 in 50 signal would occupy the frequency range 20.001 to 20.007 
| communication via satellite 410 with a hubstation 400. The GHz. However, using single-sideband modulation can make 
J hubstation is connected, for example via a local exchange, to it difficult to preserve the phase coherency between the beam 
/ a public switched telephone network (PSTN) to allow calls signals. Accordingly, double-sideband modulation tecb- 
/ to be placed between the portable phones and any telephone niques can be used instead. For example, the 1-7 MHz video 
/ subscriber worldwide, as well as between the satellite 55 signal can be used to frequency or phase modulate a 20 Ghz 
'1 phones. The satellite receives signals from the portable carrier frequency. The frequency range occupied by the 
phones at a relatively low microwave frequency, such as modulated carrier would then be approximately 
1600 MHz. At such frequencies the transmitters in battery 19.993-20.007 MHz, or more, depending on the frequency 
operated phones can be efficient and their antennas can be or phase deviation employed. To allow some margin over the 
small and omnidirectional. The satellite translates the 60 bare 14 Mhz bandwidth, a 25 Mhz carrier spacing might be 
received signals from 1600 MHz to a higher frequency for appropriate in the 20 Ghz band. Thus, 37x25 or 925 MHz 
relaying to the hubstation. can be used for the one-way satellite-hub link. This band- 
A higher frequency can be used because the bandwidth width can be halved by intelligent use of orthogonal polar- 
needed on the satellite-to-bub link is at least n times the izations as described later. ■ 
bandwidth allocated at 1600 MHz for each beam, where n is 65 FIG. 8(a) shows an exemplary satellite transponder for I 
the number of beams. For example, if 6 MHz of bandwidth the hub-mobile relay direction. The same method described j 
is re-used in each of 37 beams at 1600 MHz, then at least above for the mobile-to-hub transmissions can be used in ' 
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reverse for the coherent transport of multiple beam signals may be applied to further compensate for distortion in the 

to the satellite. The hub station (not shown) transmits a associated transponder channel power amplifiers. Moreover, 

number of Ka band frequency or phase modulated carriers to in the case of the over-dimensioned matrix power amplifier 

the satellite. These are received using a suitable Ka band described in the above-incorporated "Waste Energy Control 

antenna 360, amplified in a common low-noise amplifier 5 and Management in Power Amplifiers" application, the 

350, and then fed to FM receiver bank 340 where each number of feeder links is greater than the number of inde- 

carrier is demodulated by a respective receiver to produce a pendent beam signals to be created, thus affording a measure 

video frequency version of the signals for transmission in of redundancy against failure. 

respective beams. These video signals, for example occu- FIG. 9 shows a block diagram of a hubstation according 
pying the band 1-7 MHz, are then unconverted in respective 10 to an exemplary embodiment of the invention. The hub 
upconverters 320, using a common local oscillator 330 to antenna 600 receives Ka band carriers from the satellite and, 
preserve relative phase relationships, and then amplified after common low-noise amplification and optional down- 
using power amplifier matrix 310 for transmission via conversion in block 610, the signal is divided between a 
multi-beam antenna 300 to the mobile phones. A suitable number of receivers for respective Ka-band carriers to 
frequency for the satellite-to -mobile link is, for example, 2.5 IS obtain the beam signals B 2 . . . B„. Each beam signal is 
GHz (S-band). The amplifiers in the power amplifier matrix composed of a multiplicity of voice-modulated channel 
can be linear amplifiers to reduce intermodulation between frequencies which are separated in channel separation filters 
signals destined for different phones. The power amplifier 630. 



f matrix can for example, either be a bank of n separate 
amplifiers each associated with respective beams, or a bank 
r of N (greater or equal to n) amplifiers coupled by nxN Butler 
^matrices at their inputs and Nxn Butler matrices at their 



The channel separation filters 630 can be analog compo- 
20 nents such as crystal filters, and may involve frequency 
conversion of a selected channel frequency to a common 
lower frequency (e.g., 12.5-25 KHz, or 455 KHz) for A/D 



outputs. The effect of the Butler matrices is to use each conversion. The selected channel signal having been con- 
amplifier to amplify part of every beam signal, thus evening verted to a suitable frequency is A/D converted in A/D 
the load, providing graceful degradation in the event of 25 convertors 640. An exemplary A/D converter technique 
failure, and reducing intermodulation by absorbing a pro- suitable for use at low intermediate frequencies such as 455 
portion of the intermodulation energy in N-n dummy loads. KHz is the technique described in U.S. Pat. No. 5,048,059 
Examples of such power amplifier matrixes can be found in to Paul W. Dent entitled "Log-Polar Signal Processing", 
U.S. patent application Ser, No. 08A79,947 entitled "Waste which is incorporated here by reference, which preserves the 
Energy Control and Management in Power Amplifiers'* and 30 full complex nature of the signal by simultaneously digitiz- 
filed on Jan. 11, 1994 which is incorporated here by refer- ing its phase and its amplitude. Instantaneous phase can be 
ence. digitized for example, using the technique described in U.S. 

According to another exemplary embodiment of the Pat. No. 5,084,669 to Paul W. Dent entitled "Direct Phase 

present invention, in communication systems using TDMA Frequency Digitization", which is also incorporated here by 

signals relayed through an earth -orbiting satellite having a 35 reference. Phase digitization of all n beam signals corre- 

communications transponder using such a matrix power sponding to one channel frequency can be carried out using 

amplifier, the power amplifier can have its input Butler the technique described therein by repeating certain ele- 

combining network located at the ground station instead of ments (i.e., the trigger circuits and holding registers) n times 

the satellite. A Butler combining operation may be per- and sharing others (i.e. s the reference frequency counter) as 

formed by digital signal processing at the ground station to 40 necessary to preserve relative phase coherency. Alternately, 

form weighted sums of the desired beam signals to generate digital filters can be used instead of analog filters if the 

drive signals corresponding to each amplifier of the matrix composite beam signals are first digitized, in which case the 

power amplifier. These weighted sums are transmitted using A/D converters 640 in FIG. 9 would not be needed, 

coherent feeder links to the satellite's communications tran- The numerical results of A/D conversion are fed sample 

sponder which receives them and translates them to a second 45 by sample to numerical matrix processor 650. There is one 

frequency band for driving the power amplifier in such a such processor per channel frequency, but only the processor 

way that, after Butler combining the power amplifier for channel frequency (m) is shown for clarity. The matrix 

outputs, the output signals correspond to signals desired to processor processes the digitized beam signals to separate 

be transmitted in different antenna beam directions to out up to n separate mobile phone transmissions M 1 . , . M„ 

respective ground-terminals, which may be, for example, a 50 and transfer a sample stream corresponding to each mobile 

small handportable station. phone transmission to voice channel processor 660. The 

The resulting satellite circuitry is shown in FIG. 8(a). voice channel processor numerically performs demodulation 

Note that the input combiner which is normally present has of the signal and error correction decoding and transcoding 

been omitted since this function is now performed at the of digitized voice from the bit rate and format used over the 

ground station, as illustrated by the dashed rectangular 55 satellite to standard PCM telephone system format for 

outline 800. The antenna 810, signal processing including connection via a digital exchange (not shown) to the PSTN, 

linear amplifier 820, feeder link receivers and down convert- Thus the exemplary structure shown in FIG. 9 accomplishes 

ers 830 and output combiner 840 can be implemented in the decoding of nxm voice channels, where n is the number of 

conventional manner and thus are not further described beams and m is the number of frequencies per beam. For 

herein. 60 example, with n-37 and m-400, the system has a 14800 

This embodiment may be advantageous for certain voice channel capacity potential, 

situations, for example, dynamic reallocation of power The explanation of FIG. 9 relates to a system wherein one 

between antenna beams and timeslots may be accomplished voice channel is carried per frequency (i.e., a Frequency 

without large variations in the corresponding forward feeder Division Multiple Access (FDM A) system). However, the 

link signals, because each feeder link carries part of every 65 present invention can also be applied to Time Division 

beamsignal instead of all of one beam signal. Additionally, Multiple Access (TDMA) systems. In TDMA systems, sev- 

pre -distortion of signals sent on the forward feeder links eral mobile phone signals are carried on the same channel 
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frequency by dividing a repetitive frame period into time signal from the first antenna feed element's (beam) down- 
slots, and allocating one time slot in each frame to one convenor (not shown) is fed to the voltage control input of 
mobile phone signal. The exemplary block diagram of FIG. a 20 GHz voltage controlled oscillator (VCO) 1000. The 
9 can be applied even when the sample streams from AID video signal frequency modulates the VCO. Successive 
converters 640 represent TDMA signals. However, the 5 VCO's with their center frequencies offset by the desired 
matrix processor 650 will now separate a different set of channel spacing (e.g., 25 MHz) are used for the signals from 
mobile signals in each timeslot, so that the matrix coeffi- antenna feed elements 2,3 . . * to n/2. The VCO center 
cients are now multiplexed between several sets, each of frequency for signal 2 according to this exemplary embodi- 
which correspond to a time slot. This can be an economic ment is 25 MHz (i.e., 1 GHz/40) higher than that for signal 
arrangement because, for a given number of voice channels 10 1, and the VCO 1001 frequency for signal n/2 is thus 
complexity, the channel filters 630 will be fewer in number (n/2-l)xl GHz/40=(n-2)/80 GHz higher than that for signal i 
by a factor equal to the number of timeslots per carrier, the 1. The signals from the VCOs are summed in summer 1002 j 
A/D convenors are correspondingly fewer, the number of which can be, for example, a waveguide or stripline direc- 
matrix processors is reduced correspondingly although each tional coupler network, and the sum is amplified in a 
has to operate at a higher input sample rate, and each voice is common amplifier 1003 which can, for example, be a 
channel processor can sequentially process the signals in travelling wave tube amplifier (TWTA). 
consecutive timeslots and thus achieve the same total num- A parallel arrangement is used to deal with the other half 
ber of voice channels capacity while economically time- of the video signals numbered n/2+1 ton. The VCO 1004 for 
sharing components. signal n/2+1 is offset by half a channel spacing (i.e., by 12.5 
Each numerical matrix processor 650 is shown receiving 20 MHz in the above example where channel spacing is 25 
a control signal. This control signal can be generated by a MHz) from that of VCO 1000 and this offset is maintained 
separate computer (not shown) which controls the connect- up to VCO 1005 such that the set of frequencies used in the 
ing and disconnecting of calls to mobile phones, requiring parallel arrangement are all offset half of a channel from 
changes to the matrix of coefficients used by the processor those of the first arrangement. This minimizes any interfer- 
for separating out mobile signals from the beams. It was 25 ence which may be caused by imperfect polarization isola- 
mentioned earlier that this separation can be achieved if the tion in the dual polarization transmission system. The out- 
inverse of the C matrix was not numerically ill-conditioned, puts of the two TWTAs are connected to, for example, 
If two mobiles are located exactly at the same point on the dual-circular-polarized horn antenna 1009 via a polarizer 
earth, their two corresponding columns of the C matrix will 1008. The function of the polarizer 1008 is to launch a Right 
be identical, which causes the determinant to be zero and the 30 Hand circularly polarized signal into horn antenna 1009 
inverse not to exist. Thus, for the C-matrix to be invertible corresponding to the signal from TWTA 1003 and simulta- 
the mobiles shall be spaced far enough apart on the ground. neously a Left Hand Circularly polarized signal correspond- 
If they approach each other too closely, the C-matrix ing to the signal from TWTA 1007. 

becomes ill-conditioned. At the hubstation, the composite signal is received by a 

I According to an aspect of the present invention, however, 35 dual circularly polarized antenna and the two polarizations 

|1| when this situation threatens, one of the two (or more) are split into two respective banks of FM receivers. The 

f approaching mobiles changes frequency to a channel where center frequencies of the FM receivers correspond to the 

the other mobiles using the same frequency are adequately center frequencies of the VCOs of FIG. 10. The demodu- 

separated. It is the function of the control computer to lated outputs from the FM receivers reproduce the signals 

determine, at least at call set-up and optionally at intervals 40 received at the satellite L-band antenna elements preserving 

thereafter, which of the available channel frequencies is their phase and amplitude relationships. FIG. 11 shows an 

most suitable for allocating to a new mobile, or for handing exemplary relationship between the K-band transmission 

over an ongoing conversation. If there is no free capacity in spectra for the two polarizations, showing how the half- 

a system the system is said to be blocked and subscribers channel offset between the RHC and LHC center frequen- 

cannot place calls, much to their annoyance. When the 45 cies minimizes interaction. 

system is underloaded there are, at least on some Those skilled in the art will readily appreciate that the 

frequencies, fewer mobile signals than beams, thus the block diagrams of FIG. 10 is merely illustrative of an 

matrix C is not square. It will be shown later how the excess exemplary arrangement of coherent signal transmission 

degrees of freedom provided in underloaded systems can according to the present invention and that many functional 

then be used, not only to separate mobile signals from each 50 equivalents flow therefrom. For example, it might be advan- 

other thus avoiding mutual interference, but also to maxi- tageous first to generate the frequency modulated signals at 

mize the signal quality received from the worst-case mobile. a lower frequency than 20 GHz, for example 2-3 GHz and, 

This solution changes when an extra mobile signal has to be after summing, to convert the composite signal to 20 GHz by 

accommodated and the control computer can evaluate in mixing the summed signal with a common 18 GHz local 

advance the potential impact on signal qualities. Thus a 55 oscillator and selecting the upper sideband with a bandpass 

strategy for allocating a channel according to an exemplary filter. 

embodiment of the present invention is to evaluate the The above discussion has centered on the coherent trans- 
impact on the signal quality corresponding to the worst case portation of signals received by the satellite L-band antenna 
mobile on each channel through the inclusion of the new elements to the hubstation. The same function, namely the 
signal in the computations. The channel which suffers the 60 transportation of signals generated at the hubstation, is used 
least degradation, or meets the highest quality for the in reverse for radiation by respective satellite antenna 
worst-case mobile, is then logically selected as the one to elements, e.g., by the transponder of FIG. 8. The hubstation 
use for the new signal. This results in the group of mobiles can use a similar arrangement to FIG. 10, but with a set of 
assigned to the same frequency being as widely spatially K-band frequencies different from those used in the satellite- 
separated as possible. 65 to-hub direction, and with a larger antenna at the ground end. 

FIG. 10 shows an exemplary arrangement for coherent The satellite can employ a second, dual-polarized horn 

transmission of "video" signals from each beam. The video antenna for reception, or alternately use the same horn 
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aatenna and polarizer 1008 and 1009 as in FIG. 10 with the imaginary axes, and preserving the correct amplitude re la- 
addition of a transmit/receive duplexing filter for each tionships between the I and Q signals will then preserve the 
polarization to separate the transmit and receive signals. vector relationships including relative phase. The 2n I and Q 
Linear amplifier 350 can be duplicated for each polarization video signals can be used to frequency modulate 2n K-band 
and used to feed respective halves of FM receiver bank 340. 5 carriers having less than half the channel spacing previously 
The same half -channel frequency offset between the carriers used in FIG. 10, for example 10 MHz. While it may appear 
of the two polarizations is also advantageous in the hub-to- more spectrally efficient at K-band to use this method, it is 
satellite direction. difficult in practice to handle I,Q signal components down to 

FIG. 12 shows an alternative arrangement according to true zero frequency due to DC offsets and frequency errors, 
another exemplary embodiment of the present invention for to Consequently it is desirable to employ AC coupling and 
coherently transporting multiple signals between the hub thus exclude a portion, for example 0-50 KHz, of the 0-3 
and the satellite. In this figure, each satellite transponder MHz video signals from transmission. This places a notch in 
channel corresponding to one antenna feed element is shown the center of the 6 MHz wide L-band bandwidth that is 
as a double downcon version process comprising an antenna transponded by this exemplary method. Depending on the 
filter 1200, a low-noise amplifier 1201, image rejection filter is nature of the signals, this notch may be of no consequence. 
1202, first downconvertor 1203, IF filters 1204, 1206, IF For example, in copending, commonly assigned U.S. patent 
amplifier 1205 and quadrature downconvertors 1207, 1208. application Ser. No. 08/179,954 entitled "TDMA/FDMA/ 
The first downconvertors 1203 can use the same local CDMA Hybrid Radio Access Methods" and filed on Jan. 11, 
oscillator signal for all channels to preserve relative coher- 1994, which is incorporated here by reference, there is 
ency. The quadrature downconvertors 1207 and 1208 can 20 disclosed a hybrid access method suitable for satellite- 
use the same second local oscillator cosine and sine refer- cellular applications in which signals are conveyed on the 
ence signals for all channels, again to preserve relative downlink (satellite-mobile) by wideband TDMA in which 
coherency. The quadrature downconvertor outputs, for each mobile signal occupies an assigned timeslot in a 
example in the 0-3 MHz range, are split in cross-over repetitive frame structure, and on the uplink (mobile-to- 
network 1209 into 0-50 KHz components on lines 1215 and is satellite) by Frequency Division Multiple Access (FDMA) 
1216 and 50 KHz-3 MHz components on lines 1217 and or a combination of FDMA and Code Division multiple 
1218. The 50 KHz to 3 MHz components correspond to Access (CDMA). For example, a 63536 megabit per second 
uplink traffic channels using, for example, FDMA, an TDMA signal comprising 512 timeslots can be transmitted 
FDMA-plus-narrowband-CDMA hybrid or narrowband from the hubstation for transponding through the 6 MHz 
FDMA/TDMA, and are used to modulate separate I and Q 30 bandwidth of each satellite antenna feed element to a 
K-band transmitters for relaying these signals coherently to corresponding number of mobile phones in each cell. The 
the hubstation. These components modulate the I and Q omission of a small fraction of the bandwidth in the center 
voltage controlled oscillators 1210 and 1211. The outputs of of the channel will not disturb the character of such a signal, 
these oscillators are summed in a K-band summing network significantly, and such disturbance as may occur can be 
and the sum fed to a common TWTA for amplification to the 35 compensated at the receiving radio using the technique for 
desired downlink transmit power level. It is also advanta- DC offset compensation disclosed in commonly assigned 
geous to combine half the VCOS, e.g., the I VCOs into a first U.S. Pat. No. 5,241,702 to Paul W. Dent entitled "D.C. 
TWTA to form a signal for transmitting using RHC Offset Compensation in a Radio Receiver" which is incor- 
polarization, the other half being transmitted with LHC A porated here by reference. 

similar arrangement can be employed at the hubstation for 40 When such 512-timeslot TDMA formats are used on the 
coherently conveying the composite signal for each beam to downlink, one or more timeslots can be dedicated for use as 
the satellite. a common signalling channel, also known as a calling 

The corresponding K-band receivers would comprise an channel, forward control channel, or paging channel. The 
FM receiver for each of the I signals and an FM receiver for calling channel is used by the system to broadcast calls to 
each of the Q signals. These FM receivers would preferably 45 mobile phones originating from the network (e.g., from a 
have automatic frequency control (AFC) which removes DC PSTN subscriber or from another mobile phone). When a 
and low frequency components of the I and Q signals',. mobile detects its own phone number or ID in such a 
equivalent to having a notch in the frequency response in the - broadcast message, it replies using a corresponding uplink 
center of the channel. This is of little consequence for channel commonly called the "random access channel". The 
wideband TDMA signals and for FDMA simply means not 50 random access channel is so called because it is also used by 
using the channel in the center of the band for traffic. mobile phones to place mobile-originated calls, that is to 

In the satellite, the outputs of the K-band receivers are request service from the network. With a large population of 
reconstituted I and Q signals that are used to modulate COS roaming mobile phones, these requesting events seem to the 
and SIN L-band carriers using a quadrature modulator to system to arise more or less at random, 
produce coherent beam signals. These are applied to L-band 55 According to the aforementioned "TDMA/FDMA/ 
power amplifiers for each beam or to the PA of the afore- CDMA Hybrid Access Methods" patent application, there is 
mentioned matrix type. associated with each downlink timeslot a corresponding 

The frequency arrangements used can be similar to those uplink carrier frequency. Thus to employ the aforementioned 
depicted in FIG. 11, with RHC polarization being used, for disclosure in conjunction with the I, Q version of the present 
example, for I components and LHC polarization for Q 60 invention, the uplink carrier frequency associated with the 
components, and the carrier spacings being reduced such downlink calling channel timeslots can be chosen to corre- 
that they are commensurate with a 3 MHz modulating signal spond to the +/-50 KHz in the center of the 6 MHz 
instead of 7 MHz. A half -channel offset between the RHC- bandwidth and is used as the random access channel, 
and LHC-polarized carriers is also advantageous in this I,Q Accordingly, the 0-50 KHz signals from the crossover 
method. 65 network 1209 represent random access signals and because 

Thz I and Q signals represent, respectively, the projec- of their relatively low bandwidth the option of digitizing 
tions of the complex received signal vectors on the real and on-board and transmission by digital means to the hubstation 
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exists. This is carried out by A/D convenors 1212, the MHz of the original N signals. A suitable modulation 

outputs from each channel of which are multiplexed in scheme can, for example, be 16 QAM. In 16 QAM, four bits 

multiplexer 1213 to form a composite bitstream on the order of data are conveyed per transmitted symbol, by mapping 

of 60 MB/s which modulates a digital transmitter 1214 for two bits to one of four K-band carrier real vector values (i.e., 

transmission to the hubstation. 5 the amplitude of a cosine carrier component) and two bits to 

According to yet another exemplary embodiment, antenna one of four imaginary vector values (i.e., the amplitude of 

element signals can be transported coherently between the the sine component). The 4x4 grid of possible points that 

ground station and the satellite without bandwidth expan- result is shown in FIG. 14. Using 16 QAM, I bitpairs are 

sion. FIGS. 13 and 14 illustrate an exemplary coherent mapped to the K-band I axis and Q bitpairs to the K-band Q 

transmission method and apparatus which is based on analog 10 axis using DtoA convertors 1310 and 1311. Finally, the 

to digital conversion of each of the antenna signals followed desired K-band vector components are formed by applying 

by digital multiplexing and then modulation of the multi- the outputs of two-bit DtoA convertors 1310 and 1311 to a 

plexed stream on to the K-band feeder link carrier by means K-band Quadrature Modulator 1312 which is driven by 

of Quadrature Amplitude modulation. FIG. 15 illustrates an K-band cosine and sine carrier waves (not shown) to form a 

alternative apparatus derived from FIG. 13 which equates to 15 modulated output signal for transmission via the K-band 

infinite AtoD and DtoA precision, thus permitting the feeder link antenna (also not shown). 

AtoD's and DtoA's of the exemplary embodiment of FIG. The multiplexer can preferably have more inputs than 

13 to be replaced by analog multiplexing. signals from antenna channels. For example, a typical 

With reference to FIG. 13, operation of this coherent antenna arrangement can be a hexagonal array of 61 antenna 

transmission system is as follows. A 2 GHz signal received 20 elements. A 64-input multiplexer can then be suitable, as a 

from one of a plurality of satellite-borne antenna elements is power of two arises naturally in multiplexer construction, 

low-noise amplified and downconverted against cosine and The 3 spare inputs can then be connected to reference I,Q 

sine local oscillator signals using mixers 1301 and 1302. If signals equal respectively to (0,0), (1,0) and (0,1). The 

the bandwidth at 2 GHz that is downconverted is 5 MHz, ground station receiver can use these reference signals to 

then the resulting I and Q signals are of bandwidth 2.5 MHz 25 synchronize its demultiplexing and to determine quadrature 

each. Thus the desired bandwidth of 5 MHz may be imposed modulator carrier leakage (offset) from the (0,0) case, and to 

by the use of 2.5 MHz cut-off low pass filters 1304, 1305 provide phase references from the (1,0) and (0,1) cases for 

operating on the I and Q signals. These mixers, filters and discriminating the I-axis bits from the Q-axis bits. 

AtoD convertors 1306, 1307 are repeated for each separate In case two-bit quantization is inadequate, AtoD conver- 

antenna element signal so treated. The mixers can receive 30 tors 1306 and 1308 can be of a higher resolution, for 

the same local oscillator signals cos(wt) and sin(wT) so as example four bits. Then each 4-bit I and 4-bit Q sample will 

not to introduce any relative phase shift between channels. represent one of 256 possibilities, and this can be transmitted 

The baseband I and Q signals after filtering are converted using 256 QAM in the same way as described above for 16 

using AtoD convertors 1306 and 1307. These are arranged to QAM. However, a simplification is possible by noting that 

sample and convert the I and Q signals at least at the Nyquist 35 the complementary operations of AtoD conversion per- 

rate, which is twice the bandwidth or, in this example, 5 formed in blocks 1306, 1308 and mapping to a symbol 

MS/S. Sampling at least at the Nyquist rate allows the performed by DtoA convertors 1310 and 1311 simply cancel 

signals to be faithfully reconstructed from the samples. By each other out and can be omitted in this alternate exemplary 

way of example, the AtoD convertors are illustrated as embodiment. Then, the full unquantized accuracy of the 

having only two bits resolution, that is each I or Q signals 40 analog I and Q signals from the low pass filters is preserved 

is classified as lying nearest to one of the four values -3,-1, through the multiplexers and the digital multiplexers are 

+1 or +3 arbitrary units, as indicated by a digital code 11, 10, replaced with analog multip lexers as shown in FIG. 15. 

01 or 00. In FIG. 15, baseband signals are produced by downcon- 

In certain applications, two bits quantizing may indeed be vertors 1501, 1502 and low pass filters 1503, 1504 as 

sufficient. Such applications are characterized by the total 45 described above with respect to FIG. 13. The I,Q signals are 

signal- to -noise ratio in the 5 MHz bandwidth at 2 GHz being however no longer digitized and instead are applied directly 

very low or even negative. This can arise, for example, when to the inputs of analog multiplexers 1505, 1506 along with 

the signal bandwidth has artificially been widened by the use corresponding signals from other antenna channels (not 

of coding or spread-spectrum techniques. If the signal-to- shown). The multiplexed I samples then modulate a K-band 

noise ratio is low or negative, a few bits resolution suffice to 50 cosine carrier and the multiplexed Q samples modulate a 

make the digital quantizing noise lower than the radio noise K-band sine carrier, by use of quadrature modulator 1507. 

to avoid degradation. Those skilled in the art will appreciate Spare inputs of the analog multiplexers, as previously 

that for applications having higher signal-to-noise ratios, indicated, can be used to multiplex and transmit reference 

more bits can be used to provide greater precision. values such as (0,0), (1,0) and (0,1) which can be helpful in 

With the two-bit example, bit-pairs representing instan- 55 assisting the ground station receiver to acquire demulti- 

taneous I samples and Q samples are collected from all plexor synchronization and in correcting certain errors in the 

antenna elements and multiplexed using digital multiplexers quadrature modulator such as carrier imbalance (carrier 

1308 and 1309. The output of digital multiplexers 1308 and leakage, offset) and imperfect quadrature (i.e., the cosine 

1309 is a two-bit I and two-bit Q signal, respectively, for and sine carriers are not exactly 90 degrees apart), 
antenna number 1, followed by the same for antenna number 60 The configuration illustrated in FIG. 15 has an advantage 
2, then 3, 4, etc., until antenna 1 is again sampled. The that substantially no bandwidth expansion of the signal takes 
succession of two-bit values of I and Q is then to be place from 2 GHz to K-band. The N, 5 MHz wide antenna 
transmitted by modulation onto the K-band feeder link signals received at 2 GHz are retransmitted at K-band using 
carrier frequency. substantially the same 5N MHz bandwidth. Furthermore, no 

Since the number of bits per second is 4Nx5MS/S=20N 65 quantization noise is introduced. 

Mbits/S, a bandwidth-efficient digital modulation scheme is A suitable analog multiplexer for the exemplary embodi- 

required to avoid the signal occupying more than the 5N ment of FIG. 15 can be constructed as a binary tree, in which 
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pairs of 5 MS/S signals are first multiplexed in relatively modulating a K-band feeder link or to digitize them using 
low-speed, 2-input multiplexers to form 10 MS/S signals. the method of VS. Pat. No. 5,048,059, which was earlier 
Then pairs of these are multiplexed in higher speed 2-input incorporated by reference, prior to multiplexing, 
multiplexers to form 20 MS/S signals and so on. The FIG. 16 shows transmit signal processing in the hubsta- 
multiplexers can be constructed in a bipolar, CMOS or 5 tion for this exemplary embodiment of the present invention. 
BiCMOS integrated circuit using current steering in which Each voice channel to be transmitted to a mobile phone can 
a signal is applied to the junction of two transistor inputs be received either as a standard 64 KB/s PCM signal or as 
(e.g., emitters) that are alternately enabled or disabled by a an analog signal which is converted to PCM. The PCM 
control signal (applied, e.g., to bases) to either pass the signal is then transcoded to a lower bit rate, such as 4.8 
signal current through one of the devices or to shunt it away 10 KB/s, using a conventional voice compression algorithm 
through the other. Gallium arsenide technologies, such as such as CELP (codebook excited linear prediction), RELP, 
HBT, are also very suitable for constructing high speed VSELP or sub-band coding. The transcoded voice signal is 
multiplexers. then subject to error correction coding and supplementary 
The ground station processing system receives the time- bits can be added such as the Cyclic Redundancy Check bits 
multiplexed antenna signals on K-band, converts those 15 (CRC), Slow Associated Control Channel signalling infor- 
signals down to I,Q baseband signals of 2.5N MHz band- mation (SACCH), per-slot syncwords and inter-slot guard 
width each, and then demultiplexes them into N, separate 2.5 symbols. This per-channel processing takes place in voice 
MHz bandwidth signals of 5 MS/S each (the Nyquist rate or processing channel cards 1600. The output bitstreams from, 
higher). These signals can then be digitized on the ground to for example, 500 such channel cards, are then multiplexed 
whatever accuracy is required for further processing such as 20 with a control channel data stream from a control processor 
using an equalizer for removing inter-sample interference on (not shown) in multiplexer 1601 to form the TDMA 
a sample caused by smearing from adjacent samples due to bitstream, for example, of 6.5536 megabits per second. This 
deliberate or accidental bandwidth restrictions in the K-band is submitted to a digital modulator 1602 that numerically 
transmitter, receiver or propagation path. Such an equalizer converts the information stream to a stream of complex 
operates by subtracting a defined amount of a previous and 25 numbers at a sample rate of, for example, eight samples per 
subsequent complex (I,Q) sample value from a current bit representing the I, Q components of a modulation wave- 
value, the defined amounts being given by complex coefE- form. 

cients that are chosen to cancel inter-sample interference. The TDMA signal produced as described above is tar- 

This process can also be applied in reverse for conveying to geted for transmission to a first set of, for example, 500 

the satellite using K-band complex signal vector samples for 30 mobile phones in a particular cell or area. A number of other 

transmission by respective antennas at, for example, S-band. such TDMA signals formed by similar circuitry 1600, 1601, 

The N separate, complex (I,Q) sample streams are first 1602 are produced for transmission to other sets of 500 
preferably subjected to pre-equalizing at the ground station mobile phones in 36 other cells. The total number of cells 
such that they will be received with zero intersample inter- (e.g., 37 in this exemplary embodiment) times the number of 
ference at the satellite. Then the time -multiplex modulated 35 traffic channels per cell (e.g., 500) gives the total system 
K-band signal is downconverted in the satellite against a capacity as 18500 voice channels. The signals in timeslots 1 
K-band local oscillator to give multiplexed I and Q streams. of each cell are transmitted simultaneously on the same 
If desired, two or more stages of downconversion can be frequency to their respective cells. To avoid spillover inter- 
employed so that amplification takes place at convenient ference from adjacent cells using the same frequency at the 
intermediate frequencies. This may also apply to the 2 GHz 40 same time, this exemplary embodiment of the present inven- 
downconverters of FIG. 13, but note that the same local tion includes matrix processor 1603 to process the signals 
oscillator signals should then be employed in all correspond- from modulators 1602 by weighted addition using a matrix 
ing stages of downconversion for each antenna element so as of 37x37 complex coefficients for each timeslot. The 37x37 
not to introduce relative phase shifts. coefficients for each timeslot are contained in coefficient 

The multiplexed I,Q streams received by the satellite can 45 memory 1605 which can be distributed within the compo- 

be demultiplexed using the same multiplex clock (not nents of the numerical signal processor but which is collec- 

shown) used for multiplexers 1505 and 1506. The onus is tively identified as a separate block 1605 in FIG. 16 to better 

thus on the ground station to transmit a signal taking into illustrate the principle. During the first timeslot, a first set of 

account propagation time such that the signal will arrive in coefficients C is selected from the memory and used to 

the correct timing relationship to ensure proper demultiplex- 50 matrix -multiply the modulation signals from modulator 

ing on board the satellite. In this way, the satellite function 1602 to obtain signals for D/A convertors 1604. Each D/A 

is kept simple and reliable and complexity is restricted to the convenor can be a dual-channel unit capable of operating 

ground, where equipment can be repaired should it fail. with complex numbers. For example, the output signals - 

The above description has been simplified for purposes of from the matrix processor can consist each of a 12-bit real 

illustration to the case where all antenna element signals are 55 (I) and 12-bit imaginary (Q) part which are D/A converted 

time-multiplexed to a single TDM complex sample stream. to produce analog I, Q signals. The I, Q signals are fed to FM 

Those skilled in the art will readily appreciate, however, that K-band FM transmitters for transmission from the hubsta- 

a hybrid TDM/FDM scheme could be used in which groups tion to the satellite. 

of time-multiplexed signals are formed and used to modulate When transponded by the satellite to the ground on 

separate FDM carriers. This modification could be used if, 60 S-band, the result of the matrix processing will be that each 

for example, a single multiplex stream would result in an mobile phone receives only its own signal, the inter-cell 

unpractically high sample rate. interference from other cells having been cancelled by the 

It is also for the purposes of illustration that the above addition in the matrix processor of compensating amounts of 

description has concentrated on the Cartesian (I,Q) repre- opposite sign as determined by the coefficients retrieved 

sentation of complex signals. It is equally possible to form 65 from memory 1605. This is possible if the 37 mobiles using 

polar or logpolar representations of complex signals, to timeslot 1 in their respective cells are spatially separated, 

multiplex these signals using analog multiplexers prior to i.e., not both at the same location on the edge of their 
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respective cells. This condition can be maintained by the 
exemplary timeslot assignment algorithm feature of the 
present invention, which also provides a general channel 
assignment algorithm, and is based on maximizing the 
signal quality provided to the worst case mobile. 

A timeslot duration is typically about 40 if a 20 mS 
TDMA frame period is used. One timeslot corresponds to 
256 bit periods at 6.5536 MB/s and 2048 complex numbers 
are produced by modulator 1602 for every timeslot. After 
matrix processor 1603 has processed 2048 sets of 37 com- 
plex number inputs using the set of coefficients for the first 
timeslot, the coefficients are changed for the second and for 
subsequent timeslots to effect correct interference cancella- 
tion between corresponding sets of 37 mobiles using 
timeslots 2, 3 etc. 

If two mobiles receiving the same timeslot in different 
cells approach each other too closely during the progress of 
a conversation, the control processor (not shown) will note 
a difficulty in arriving at a suitable set of coefficients for 
interference cancellation. This is highly unlikely given the 
limited speed of landmobile phones in relation to typical cell 
size, but if it occurs, the control processor evaluates whether 
a timeslot change would be appropriate for one of the 
mobiles. The aim is to connect the mobile using a timeslot 
that no other mobile in close proximity is using. If necessary, 
a mobile even occupying an ideal (e.g., low interference) 
timeslot could be shifted to a just adequate (e.g., barely 
tolerable interference) timeslot to release its original 
timeslot to solve a proximity problem at hand. It is probably 
not necessary in practice to consider such a situation because 
with, for example, 500 timeslots to choose from, it would 
usually be possible to find a better timeslot than the timeslot 
currently threatening to cause bad signal quality. Allowing 
one timeslot change per cell per 10 seconds, for example, 
would be expected to achieve adequately rapid optimization 
of timeslot assignments and adequate adaptation to mobile 
movement. 

In fact a more rapid rate of adaptation is provided to 
handle the rate at which new calls are placed and old calls 
cleared down. With a capacity of 37 mobiles per timeslot 
and an average call duration of 3 minutes, a particular 
timeslot is vacated in some cell approximately every 5 
seconds and a new call is then assigned to that timeslot. 
Overall, given, in this example, 500 timeslots and 37 cells, 
100 timeslots spread over all the cells are vacated every 
second and reassigned. 

Such a communications system should be designed to not 
be loaded up to 100% of system capacity or the next call 
attempt will be blocked. With 500 timeslots per cell 
available, an average loading of 474 timeslots can be 
reached for a blocking probability of 1%. Thus, on average, 
26 out of 500 timeslots are unused on each of the 37 
multiplexers 1302 in this exemplary embodiment. It should 
be noted that it is immaterial which multiplexer is used to 
transmit a particular timeslot to a mobile. Whichever 
timeslot is selected, it is the choice of an associated column 
of matrix coefficients that determines that mobiles using the 
same timeslot are non-interfering. Thus if the same timeslot, 
for example number 371, is vacant on two or more multi- 
plexers 1601, it is immaterial which one is used to connect 
a new call. 

Thus the assignment algorithm executed in the control 
processor first determines which timeslot is vacant on the 
greatest number of multiplexers. This is the timeslot on 
which there are currently the least number of mobile con- 
versations. Using information from the random access 
receiver on the relationships between signals received from 
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the new mobile (i.e., the C-matrix coefficients determined by 
correlation of the new mobile's random access signal with 
all antenna element signals), the control processor evaluates 
the change needed to the set of coefficients in coefficient 
memory 1605 associated with the vacant timeslot to main- 
tain non-interference if that timeslot were to be used for the 
new signal. 

The general principles that explain how the choice of 
coefficients in coefficient memory 1605 is arrived at for an 
exemplary embodiment will now be outlined. 

As discussed earlier, for receiving signals from mobiles, 
antenna element 1 received an amount Cll of mobile signal 
Ml plus an amount CI 2 of mobile signal M2 and so on. To 
state this more generally, antenna element k received an 
amount Cki of mobile i*s signal. Assuming reciprocity, a 
signal Tk transmitted from antenna element k would be 
received in an amount Cki.Tk at mobile i, because the path 
from element k to mobile i is assumed to have the same 
attenuation and phase shift in each direction, given by the 
complex number Cki. 

Therefore the signals R received at the mobiles are related 
to the signals transmitted by the antenna elements by the 
matrix equations: 

R-C.T; 

where the superscript t indicates a transpose matrix. 

The transpose of C is used because the first index k of Cki 
multiplies the corresponding index of the T-element, while, 
in the mobile-to-satellite direction where the signals 
received at element k from mobile i are given by Cki.Mi, it 
is the second index i of C which corresponds to the index i 
of the mobile signal Mi that it multiplies. Thus the indices 
of the matrix coefficients are transposed in the satellite-to- 
mobile direction as compared to the mobile-to-satellite. 

In order to achieve non-interference, the set of signals 
transmitted from the satellite antenna elements should be 
given by: 

The inverse of the transpose is just the transpose of the 
inverse, therefore the set of coefficients contained in coef- 
ficient memory 1605 for downlink timesIot(j) are just the 
transpose of the set of coefficients associated with uplink- 
frequencyfl) in numerical processor 650 of FIG. 9, at least 
under the assumption of reciprocity. 

Reciprocity applies when the uplink and downlink fre- 
quencies are the same. Relative amplitude reciprocity 
applies if the antenna element patterns are the same on both 
uplink and downlink frequencies. Phase reciprocity does not 
apply, because relative phase depends on the small differ- 
ences in relative distance travelled by the signals to/from 
each element, divided by the wavelength and multiplied by 
360 degrees. If the wavelength is different on the uplinks and 
downlinks, then the phase relationships will be different. 
However, relative time delay differences are frequency 
independent and therefore reciprocal. Accordingly, a set of 
relative phase differences at one frequency can be translated 
to a set of time differences using a first wavelength, and then 
reconverted to a set of phase differences using another 
wavelength in order to derive a set of coefficients valid at a 
second frequency from a set known at a first frequency. 

Based on the foregoing discussion, the coefficients for 
transmit contained hi memory 1605 according to an exem- 
plary embodiment of the present invention can be deter- 
mined by the following steps: 

(1) correlating the signal received from a new mobile 
during its random access transmission with the indi- 
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vidual antenna beam element signals to determine a Substituting for T from (2) into (1) we get: 
new column of coefficients for the receive C-matrix; 
(2) determining a new inverse C-matrix for receiving 
traffic from the new mobile based on the old inverse 
C-matrix and the new column; 5 



(3) transforming the new receive C-matrix column to a , . „ , . r ^, A 

new transmit C-matrix row by scaling relative coeffi- showm f that only if C A is the MxM unit 

- - matrix I 



cient phase angles using the ratio of up- to down-link 
frequencies; and 



(4) determining a new transmit inverse C-matrix based on C A ° l ( 3 ) 
the old transmit inverse C-matrix and the new row. 

An exemplary detailed mathematical procedure which is a necessary condition. Since C is not square, we cannot 

can be used to carry out the above exemplary embodiment simply invert it and write: 

is now developed for the underloaded case, i.e., the case 15 

when there are fewer currently active mobile signals than the A=c -i 
number of antenna feed elements available on the satellite to 

communicate with them. Such spare capacity is typically w ^ * » ■ A r w w . t % T 

designed for in radio telecommunications to provide a 98% Moreover, GA.-I is a art of MxM equations that the NxM 

probability of having a free channel to serve a new call so , ft u ^ owns A must MfiU so that the MxM terms of the 

that customers are not overly irritated at call blockages. 20 P roduct mdeed ^ vc mc MxM umt matnx L 

The active mobiles are designated 1 . . . m and the signals Since N>M, the number of unknowns is greater than the 

intended to be received by these mobiles are designated Rl number of equations, so there is no unique solution to 

. . . Rm for this example. The antenna element/transponder equation (3), but a continuum of solutions. Other conditions 

channels available for communications therebetween are 25 musl be imposed to define a particular solution of interest, 

designated 1 ... n and the signals fed to the antenna ^ condition imposed here is that the total power inherent 

elements for transmission by each respective element are ^ ^8^^ the vect or of signals T is minimized, 
designated Tl , . . Tn. As before, the matrix C, this time an 

mxn non-square matrix, determines how much of each It can be verified that a particular solution of equation (3) 

transmitted signal Tk reaches each mobile as Ri, the matrix ^ A-C*(CC) _1 -U, where * signifies conjugate transpose, 

is given by the equations: 3 This can be verified by substituting the particular solution U 

for A in equation (3), obtaining: 

m - en n + cnn.. . + a« r« 

Kl - C21 • 71 + C22 • 72 . . . + C2/J • 7* CU^CCiCCy^CC^iCC)' 1 , 

Rm - Cml ■ 71 + Cm! • Tl . . . + Cmn • Tn 35 

which is clearly equal to I as required, 

or simply R =C.T in matrix/vector notation. \ general solution can be formed by adding an arbitrary 

Because C is no longer square, it has no direct inverse, so matr ix V to the particular solution found above, obtaining: 
there is no unique solution for T given by: 

40 

T-CT\R A=C(CCT l + V f 

Instead there are a continuum of solutions, as we have more 

degrees of freedom to choose T values than conditions to but this must still fulfill equation (3). 

satisfy (i.e., n>m). Substituting this value of A into equation (3) we get: 

However by imposing the condition that the mean square 45 

power fed to the antenna elements in order to create the c(C(cc)- l + v) « / (4) 

desired mobile received signals R shall be minimized, the Le ^ cc{cc) + cv - I 

unique solution below is obtained: Le f + cv - i 

r=c.(cco _1 50 Le -> cv - 0 

This equation can be derived as follows. Let R Asimrf be the Thus, V may be arbitrary only so long as it fulfills 

M-element vector of signals we wish to be received at the equation (4). It is possible for a non-zero V matrix to give 

receiving stations, and T be the N-element vector of signals identically zero when premultiplied by C as long as all V's 

applied to the transmitting antennas, where N>M. C is an M columns are orthogonal to all C's rows. The rows of C are 

by N matrix of coefficients Cik that describes how the signal 55 N-element vectors, but there are only M of them, therefore 

from transmitter antenna j propagates to receiving station i. they do not totally span their N-dimensional space. There are 

Denoting by R achieved the M-element vector of signals N-M other dimensions in that space that the rows of C do not 

actually received, we thus have project into, and the columns of V may thus consist of any 

vectors that are confined to that N-M dimensional sub-space 

^«*io«r^ r (1) 60 mat jo not project into C's M-dimensional sub-space. 

We wish to find what T should be as a linear function of Thus, the general solution of equation (3) is 
the signals desired to be received, so that the minimum total 
transmit power is consumed in the process. The linear 
combinations formed by the coefficients of an M by N 

matrix A to be found are. where U is the particular solution identified above and V 

T^Ada.rrd (2) must satisfy C.V=0 



A-U+V; 



65 
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The transmitted signals T are given by by removing the worst case mobile from the current group 

the mobile is associated with and associating that mobile 
n = ah JU+An <R2 + ...Aim Rm with a different group. This is done by recomputing the 

tz - A21 Rl +A22- R2 + . ..Aim Rm above expressions with C diminished by the row corre- 

5 spending to the worst case mobile, thus determining the 
satellite power saving that would be saved in supporting 
only the remainder in the most efficient manner. Then the 
Tn - Artl ■ Ri +An2 Rl + . . .Anm • Rm removed row of C is used to augment in turn each of the C 

matrices associated with other groups of mobiles using 
where Rl, R2, etc. are the elements of R desired . If Rl, R2, 10 different frequency channels (FDMA) or multi-carrier 
etc. are all independent signals intended for different receiv- (CDMA) or timeslots (TDMA) and the above expressions 
ing stations, there is no correlation between them so they add computed to determine the increase in power that would be 
rms-wise in the linear summing process that forms the necessary to support that mobile as a member of each of the 
T-elements. other groups in turn. If the increase in power in one of these 

Thus, the mean square value of Tl is just [A11.R1| 2 + 15 cascs * less man mc power saved by removing the mobile 
|A12.R2| 2 . . . +|Alm.Rm| 2 fr° m original group, then a frequency or timeslot han- 

likewise, the mean square value of T2 is |A21iU| 2 + dover to the new group can be performed in order to improve 
|A22.R2| 2 . . . +|A2m.Rm| 2 satellite power utilization. This procedure can likewise be 

Adding these expressions down columns that contain the used for determining which of a number of existing groups 
same Ri, we get: 20 a new mobile call should be associated with, i.e., to find the 

group that would result in minimum increase of satellite 
power -Sigma, flRip-siGMAjAijl 2 ! power used when a new call is connected. 

FIG. 17 shows an exemplary embodiment illustrating the 
Now SIGMAlAijl =SIGMA l <Aij.Aij)=SIGMA f (A'ji.Aij), interconnections between the transmit and receive matrix 
where A'ji refers to element ji in the conjugate transpose of 25 processore aa d the control processor at the ground station to 

effect the above-described interference cancelling and opti- 
Butthis value SIGMAissimply the jj diagonal term of the mum channel allocation behavior, 
whole matrix product Xjk=SIGMA i (A'ji Aik), which is the receive matrix processor 1700 receives digitized 

equation for matrix multiplying A' and A, i.e., X=A'A. signal samples f rom t he ground station RF section. The 

Now substitute A-U+V; then: 30 rece { V e processing can be structured, for example, according 

X=(U*+V).(U+V)=U'U+VV+U , V+V'U and U'V+VU=2 to the exemplary FDMA embodiment of the invention of 
Re (U'V). FIG. 9, or according to the exemplary TDMA embodiment 

Substituting U=C'(CC)-\ i.e., U'=(CCT l C, into the of FIG. 16. Moreover, an exemplary CDMA embodiment 
foregoing equation, we get U*V=(CC , )" 1 .CV=0 because can be constructed by, for example, increasing the band- 
CV«0. 35 width of the channel splitting filters and including a CDMA 

Therefore 2Re(U'V)=0 and UV+VU«0. version of the per-channel processing in the circuitry of FIG. 

Hence SIGMA,- |Aij| 2 -SIGMA,- (tUijp+|Vij| 2 ), leading to: 9. Further, exemplary embodiments of the present invention 
POWER-SIGMA 7 ORif. SIGMA I -|Uij| 2 ]+SIGMA 7 - [|Ri| 2 . ma y be constructed using the novel subtractive CDMA 
SIGMAjVijl 2 ] system described in U.S. Pat. No. 5,151,919 to Paul W. Dent 

Since the two terms involving respectively U and V can 40 entitled "CDMA Subtractive Demodulation", which is 
only be positive, the power is minimized when the choice of incorporated here by reference. These features of the present 
the arbitrary matrix V in the second term is zero. Hence, the invention also lend themselves to implementation in land- 
solution for the transmit signals that create the desired based cellular systems. 

received signals is: The receive matrix processor 1700 separates the indi- 

45 vidual channel signals by applying inverse C-matrix coef- 
T-AJij^ when: A=C{CCy ficients supplied by control processor 1702 as described 

This solution also holds for the case where N=M, for then C abovc » to eliminate or suppress co-channel interference, 
is square and the above reduces to: coefficients can, for example, be determined as fol- 

lows. 

A-cr 1 50 When M spatially separated antenna/receiver chanoels 

receive different combinations Ri of M signals Si, given by 

Applying the foregoing principles, the spare degrees of 

freedom are used not just to create co-channel interference Ri=siGMAj(Cij\Sf) (5) 

free signals at every mobile, but also to maximize the 

wanted signal values for a given total radiated power. The 5S or in matrix notation, R=C.S, then the separation of the M 
total mean square radiated power is in fact the sum of the signals has a straightforward solution 



square magnitudes of the coefficients of the matrix A defined 
by: 



R^C.S' 1 (6) 



A=C.'(C.C)- 



When the number of antenna/receiver channels N is 
60 greater than the number of signals M they receive, the matrix 
j The sum of the squares down a column of A gives the C is not square and cannot be inverted. There are a coo- 
radiated power used in communicating with a corresponding tinuum of solutions possible using any subset M of the N 
mobile. The worst case mobile, i.e., that using the most channels, but there can also be a desired unique solution, 
satellite power, can thus be identified. According to an The reciprocal problem for transmitting M signals using 
aspect of the present invention, the control processor at the 65 N transmitter channels was solved above by imposing the 
hubstation periodically examines whether total satellite additional desire to minimize total transmit power. In the 
power can be minimized (or utilization of power optimized) receiving case, we can find the desired unique solution by 
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imposing the condition of maximizing the signal to noise 
ratio. To do this, a finite amount of noise must be assumed 
to exist in the receivers. 

Before this solution is described, another solution will be 
described for solving the equations: 



01 

C21 



51 + 02S2. ..+ C1M SM 
SI +C22-S2. . . + C2M SM 



CNi ■ SI + Cm • 52 . . . + CNM ■ SM 



« R2 



RN 



(7) 



When N>M there is an excess of equations over 
unknowns. They should all be consistent and solving any 
subset M of N should yield the same answer. Due to receiver 
noise, however, which causes uncorrelated errors in the 
received values R, the equations will not all be exactly 
consistent. 

A known solution to this is the so-called least-squares 
solution. The least squares method seeks the solution which 
minimizes the RMS sum of the noise errors needed to be 
added to the R-values to make the equations consistent. 

An error vector E may be defined as 

E-C.S-R (8) 

The sum square error is then 

E'E-(CS-Ry.(C.S-R)-S'.C\C.S-R'.C.S-S , .C , Ji+R'Ji (9) 

Differentiating this expression with respect to each R value 
to obtain the gradient yields: 



g?ad(E'E)-2C'.C.S-2C , Ji 

E'E is a global maximum where grad(E'E)=0, 



(10) 



i.e, 



C'CS-CR, or S-(C'C)-\C'Ji 



(11) 



This may be compared with the least-power transmit solu- 
tion where 

A=c(ccy l 

The least squares solution for reception given above is not 
necessarily that which maximizes the quality of each signal. 
To find the solution that maximizes each signal quality we 
in turn find the best A matrix row that yields that signal. 

Separated signal Si is given by row i of A, henceforth 
written Ai, multiplied by the vector of receive channel 
outputs R, i.e., Si=Ai.R. R is given by C.S+ Noise where 
"Noise" is a vector of uncorrelated noise having components 
N1,N2 ... in the receiver channels. 



Thus 

Si^Ai.C.S+Ai.Noisc 



(13) 



The amount of wanted component Si that appears in Si is 
given by 

(Ail.Cli+Ai2.C2i+Ai3.C3i . . . +AiN.CNi).Si=Ai.Ci 
where Ci means the ith column of C. Assuming all Si 



15 



The least squares solution for the M signals is thus S-A.R 
where 



(12) 
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are transmitted with unit power, the power in the 
extracted wanted component is 



iHAi.af-AiCi.Ci'.Ai' 



(14) 



There are also, however, unwanted components in the 
extracted signal due to the other signals Sk. The sum of the 
unwanted powers for all k not equal to i is given by 



10 



h-A L CdirrLCdimAi ' 

where Cdim means the matrix C with column 
In addition, there is a noise power given by 



(15) 



i removed. 



|AilJVl| 2 +|Ai2^2| 2 . . . -AiAi'jiI 



(16) 



where n is the mean square value of each of the noise signals 
N1,N2, etc. The signal-to-noise-plus interference ratio is 
then given by 



20 



AjCiCQA 
+ Cdim * Cd 



(17) 



Mathematicians will recognize this expression as the ratio of 
Hermitian forms. The maxima and minima of such expres- 
sions as 



xux 
xvx 



are given by the eigenvalues q of V. -1 U, i.e., by the solution 
30 of det(V" 1 .U-qI)=0. The values of X which give these 

extrcma are the corresponding eigenvectors. V. _1 U in our 

case is (Cdim.C , dim+nI)"'.Ci.Ci*and X is A'. 

We now use the theorem that the eigenvalues of the 

product of an n by m matrix with an m by n matrix where 
35 n>m are equal to the eigenvalues of the product taken in 

reverse order, plus n-m zero eigenvalues. 

Using as the two matrices in question the N by 1 matrix 

(Cdim.Cdim'+nl) -1 .^ on the one hand and the 1 by N 

matrix Ci'on the other hand, the eigenvalues we need must 
40 be those of the inverse product 



(18) 



Cf (Cdim - Cdim +nl)~ l ' Ci 
1*N WxN Nxl 



45 This, however, has dimension lxl, i.e., it is a scalar, so it has 
only one non-zero eigenvalue. 



Hence 

50 ^Cr^Cdim.Cdim+rt/)" 1 .^ 



(19) 



The associated eigenvector Ai' is the solution V of an 
equation of the form 

Matrix. V= V.Eigenvalue 
(Cdim.C 'dim+nO" Kd.C i '. V= Vq (20) 

Substituting for q from equation (19), 

(Cdim.Cdim+nO-^Ct.Cr. V^VCr.iCdim.Ctim+nO^.Ci (21) 

60 

It may be verified that letting V^Cdim.C'dim+nl^.Ci 
makes the right hand and left hand sides of equation (21) 
identical. Thus, this eigenvector is the optimum solution for 
the row of coefficients Ai that extract Si from R with best 
65 signal-to-noise+interference ratio. 

If instead we set out to maximize signal to (signal+noise-t- 
interference) ratio we would get 
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At^C.C'+n/y'KCi or Ai=CC .(CC'+nf)- 1 (22) 

i.e., the whole C-matrix is used in the inversion and not 
Cdim with one column removed. The value that maximizes 
S(S+N+I) should, however, be the same as that which 
maximizes (S/(N+I) as their reciprocals differ by the con- 
stant 1 only. 

It can be shown that this solution only differs by a scalar 
factor 1/(1 +q) from the solution which maximizes S/(N+I), 
and since a fixed scaling does not change signal to noise 
ratios, it is effectively the same solution. If such Ai's are now 
derived for all i and laid one under one another to form an 
M by N matrix A, the rows Cr", being the original columns 
Ci, also lie under one another to form the matrix C\ 



A-cico+n/)- 1 

This is similar to the solution for minimum transmit power 
derived above, except that the "C" matrix here is the 20 
transpose of the transmit matrix and is M by N instead of N 
by M. That means that NxN matrix CC has rank of only 
M<N and it has no direct inverse, being singular. However, 
the addition of the noise down the diagonal through the term 
nl is the catalyst that makes the matrix to be inverted 25 
non-singular and the above solution computable. 

The solution in the transmit case provided a way to test 
how much the total transmit power, having been 
reoptimized, would have to increase to support one extra 
signal. Reciprocally, in the receive case, it can be tested how 30 
the addition of a new signal to those already received would 
affect the signal to noise ratio after re-optimization of the 
above coefficients with an extra column added to the C 
matrix. The extra column of C's in question represents the 
relative strengths and phases which the new signal is 35 
received by the N receiver/antenna channels. This is deter- 
mined while the new signal is appearing on the random 
access channel and not in conflict with other signals. 
Furthermore, random access can be made with higher power 
or more coding than for normal traffic so as to facilitate 40 
detection and decoding, p The signal is decoded and retro- 
spectively the decoded signal can be correlated with signal 
samples recorded from each of the N channels to determine 
the new C-matrix coefficients. A test is then made by 
appending the new C column to each of a number of 45 
candidate C matrices in turn associated with different groups 
of ongoing signals in order to determine the group that 
would have its worst case SNR degraded the least by 
inclusion of the new signal. This then determines the allo- 
cation of a channel to the new signal for traffic, and explains 50 
how the C matrix coefficients are arrived at a row-at-a-time 
during the random access and channel allocation process. 

The separated channel signals are processed in separate 
channel processors 1701. The channel processors can either 
be engaged to process subscriber traffic, after a call has been 55 
connected, or can be employed to search for random access 
signals from a given direction. The latter is done by com- 
bining the received signals from the satellite to form beams 
covering fixed regions of the earth from which random 
access signals may be received. The coefficients used may be 60 
chosen by the control processor 1702 to provide cancellation 
or reduction of interference from other signals on the same 
frequency from other regions so as to maximize the prob- 
ability of intercepting a random access message. The ran- 
dom access message can also be provided with an additional 65 
degree of error correction coding to maximize reception 
probability in the absence of a-priori knowledge of the 
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direction from which an access attempt is received. 
Optionally, the random access channel can be frequency- 
planned to avoid immediate frequency re-use in adjacent 
cells, for example by the use of a 3-cell frequency or timeslot 
re-use plan, since using three frequencies or limes lots for 
random access does not have such a deleterious effect on 
total system capacity as if such a frequency usage plan were 
adopted for every traffic channel. 

The channel processors 1701 provide information to the 
control processor 1702 regarding the amount of each signal 
in each beam channel or separated channel, which the 
control processor 1702 uses to control the interference 
cancellation coefficients used by the receive matrix proces- 
sor 1700. Depending on, for example, the determination 
respectively of correlations between each separated signal 
and each beam signal or the determination of correlations 
between separated signals, two different control concepts 
can be applied by control processor 1702. 

In a first exemplary control implementation, a separated 
channel signal decoded by a channel processor 1701 is 
correlated or partially correlated with each non-separated 
beam signal in turn. The electrical connections for achieving 
this correlation are disposed between every channel proces- 
sor 1701 and every other channel processor 1701, however 
these connections are omitted from FIG. 17 for clarity. The 
part of the separated signal that is used for correlation can 
suitably be a known bit pattern in the channel signals, for 
example a synchronization word or bit pattern. The corre- 
lation results directly represent the C-matrix coefficients and 
these are processed by the control processor to obtain 
A-matrix coefficients as defined above. 

In a second exemplary control implementation, a sepa- 
rated channel signal decoded by a channel processor 1701 is 
correlated with at least part of other channel signals to 
determine the residual amount of non-cancelled interference 
present due to other channel signals. That part of the other 
channel signals with which correlation is performed can 
suitably be a known pattern contained in each signal, such 
as a synchronization word. Since these patterns are known 
it is not necessary to cross-couple the channel processors 
1701 to each other, thus avoiding a mass of interconnections. 
Furthermore, since adaptation of the receive matrix coeffi- 
cients by the control processor 1702 does not have to take 
place at a rapid rate, as they are relatively static for a given 
set of transmitting mobile phones, the correlation with 
different signals can occur at different times at which the 
transmitters, by prearrangement, insert a special sync word 
for the purpose of correlation. 

For example, suppose a known, 16-bit sync pattern is 
employed within each segment of transmitted signal, e.g., a 
TDMA timeslot. There are 16 possible orthogonal 16-bit 
words, so 16 different signals can be allocated orthogonal 
sync words. A Fast Walsh Transformer such as the one 
described in U.S Pat. No. 5,357,454 entitled "Fast Walsh 
Transform Processor", which is incorporated here by 
reference, provides an efficient means to correlate a signal 
simultaneously with all possible orthogonal codewords and 
thus directly determine the residual, non-cancelled interfer- 
ence amounts. If however the number of signals whose 
residual interference contributions are to be discriminated is 
greater than 16, for example 37, then 15 at a time can be 
arranged to use different orthogonal codewords while the 
other 22 use the 16th codeword. The 15 which are chosen to 
use different codewords can be changed between successive 
TDMA frames such that after slightly over two frames all 
signals have been uniquely discriminated. 

This exemplary procedure can also be applied to FDMA 
or CDMA uplink modulations. In the CDMA case, for 
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example, orthogonal spreading codes can be allocated to 
facilitate discrimination. If a hybrid FDMA/CDMAuplink is 
used with, for example, four overlapping, orthogonal 
CDMA signals on each frequency channel as described in 
the aforementioned disclosure, U.S. patent application Ser. 
No. 08/179,954, entitled "TOMA/FDMA/CDMA Hybrid 
Radio Access Methods" and filed on Jan. 11, 1994, then the 
system can readily search simultaneously for known sync 
patterns employing all four orthogonal codes. By permuting 
the underlying sync patterns as described above, it is pos- 
sible to discriminate residual interference contributions from 
any number of different CDMA transmissions using the 
same channel frequencies at different locations. This can be 
accomplished, for example, after separating the signals 
using the C-matrix, a signal can be correlated with its own 
known bit pattern and the known patterns of other signals 
that should have been cancelled; results of the latter corre- 
lations yield the amount of residual, uncancelled signal and 
can be used to update the C-matrix. 

In this second exemplary implementation, C-matrix coef- 
ficients are not directly determined, but rather the residual 
interference amounts are related to errors in the A- and 
C-matrix coefficients. This relationship can be demonstrated 
as follows. 

The satellite or base station broadcasts N combinations of 
M desired signals from N transmitter/antennas. The N 
combinations should be chosen such that each of the receiv- 
ing stations receives only its intended signal, and the other 
M-l at that receiver are cancelled. The N linear combina- 
tions are preferably those derived as set forth above, which 
result in each receiving station receiving its intended signal 
only, and with minimum total transmitter power. 



73V 
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matrix A, i.e., T=ARd. 

A is in turn shown above to be preferably equal to 
C'(CC)" 1 where Cij is the propagation from transmitter/ 
antenna j to receiver i. Estimates of Cij are made at call 
set-up time for the receive direction and transformed to 
estimates for the transmit direction as described above. 
There will, however, be errors in the estimates of the Cij for 
the transmit direction that are used to compute the matrix A. 
Let us assume that the estimated transmit matrix C is equal 
to the true matrix Co plus an error matrix dC, i.e., 

C=Co+dC or Co=C-dC 

The signals Ra actually received by the receiving stations 
are given by the true C-matrix Co times the transmitted 
signals, i.e., 

Ra=Co. T=CoAJid=(C-dC)C{CC)- l Jld=Rd-dcC.C{CC)- l Rd~ 
Rd-dCAJld 

The errors dR in the received signals dR=Rd-Ra are thus 
given by 



50 
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30 

(23) 

of the error vector dR contains a part 
^ of each of the other unintended signals j. 

If the M signals contain known signals, patterns or 
syncwords, by correlating with these at a mobile receiving 
signal i, it is possible to determine the residual unwanted 
amount of signal j, and thus determine e 0 . 

The syncwords can be orthogonal so that correlation with 
all of them can be performed at the same time by means of 
an orthogonal transform such as the Walsh-Hadamard trans- 
form. If the number of orthogonal codewords available is 
less than the number of signals M, the orthogonal codewords 
can be assigned to groups of immediately surrounding 
beams or cells whose signals are most likely to interfere due 
to imperfect cancellation. A limit set of orthogonal code- 
words can be permuted between the M signals to allow 
different subsets to be resolved at a time, and all M to be 
resolved sequentially. In this way, by correlating the 
received signals Ra over the portion containing the known 
signal pattern with all orthogonal codewords the amount of 
own codeword is obtained as well as the amount of 
unwanted codewords. The amount of other codewords is 
scaled by dividing by the complex amount of own codeword 
correlation to yield the normalized error residuals e,-, that 
may then be complex-value averaged over several measure- 
ment intervals before being reported by the receiving sta- 
tions back to the transmitting stations on a reverse Slow 
Associated Control Channel. To reduce the volume of 
reporting, each mobile can at each interval restrict itself to 
reporting only the largest error its correlator determines. The 
transmitting station can optionally either assume that the 
other errors are zero at that station, or that they are as 
previously reported if no action to correct them has been 
taken in the meantime. 

The matrix E-e^- may thus be equated to the matrix dC.A 
in (23), so we have dC.A-E or A'.dC'-E'. 

This is an insufficient set of equations for the unknowns 
dC, but a unique solution exists for which the sum of the 
squares of the dC's is least, that solution being 

Moreover, if A-C^C)" 3 then ACA'A^^C, therefore dC= 
C.E'or 

dC-E.C 

Thus, given the original estimate of C and the residual 
correlation measurements reported by receiving stations, the 
error dC in the original estimate may be calculated and the 
estimate of C gradually refined. 

As mentioned above, if the reverse SACCH signalling 
capacity does not allow all errors to be reported every time, 
it is sufficient to report only the largest. The transmitter can 
choose only to correct the largest there and then, or to wait 
until others are reported. In order to ensure that others are 
reported, the transmitter can request the receiver to make 
specific measurements via the forward SACCH channel. 
These refinements are mentioned for the sake of complete- 
ness in describing the scope of the invention, but the extra 
complexity is probably not needed in a satellite-mobile 
communications system where the relative positions of 
mobiles in the satellite beams changes only slowly relative 
to the speed of communications. 

The control processor obtains initial estimates of the 
downlink C and A matrix coefficients, measured on the 
uplink by syncword correlation as described earlier, to the 
downlink frequency. The control processor then continually 
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outputs corrected A-matrix coefficients suitably translated to 7-site, 3-sector pattern employed conventionally, showing 

the downlink frequency as described above to transmit that what has been gained from using the same frequency in 

matrix processor 1704. every sector has been lost by the need to increase the re-use 

A complication can arise in performing this translation pattern size from 7 sites to 21 sites. Thus, according to this 

due to phase mismatches between each antenna element 5 exemplary embodiment of the present invention three or 

channel. It was stated above that the relative amplitude more sectors or antenna elements around the 360 degrees of 

between signals on the uplink and downlink frequencies azimuth should be used. 

could reasonably be considered to be the same, and that the FIG. 18 shows an exemplary cylindrical array of slot 

relative phase between signals can be scaled by the ratio of antennas 1800 suitable for implementing the present inven- 

up- and downlink wavelengths. However, consider the case 1Q tion in landbased cellular systems. The array consists of 

where phase mismatches exist between the channels that rings of eight slots around a metallic cylinder. Horizontal 

relay the mobile-satellite uplink signals from each antenna slot antennas give the desired vertical polarization, and the 

element. The signal phases are then not just antenna element slots m a nalf wavelength long, e.g., approximately 16 

phases, PHI®, but contain the additive mismatch terms, centimeters for the 900 MHz band. It can be desirable to 

dPHI(i). If PHI(i)+dPHI(i) is then scaled by the ratio of 1S cm Pj ov alternatively circular polarization jit the base station 

wavelengths, the PHI(i) part will scale correctly but the 15 ^^ined with hnear polarization at the mobile phone, 

mismatch part dPHI(i) will oot because there is no correla- especially when the mobile is a hand portable of uncertain 

. w . tl _ 4l _ j. .. . antenna onentation. Circular polarization can be formed by 

Hon be^een pba* mismatches on the up and downlink ^ sl cmssed £ ks Qr { e 

paths.Ifmeup-anddo^knk^^^ combination for the array elements. It is often convenient 

respectively by uPHI(i) and dPHI(i) then we need to calcu- 2Q when ^ ^ structures t0 form both ^1^^ 

: simultaneously, and this can be exploited by using opposite 

a.(Pffl(*)-«PHi(0>HiPHi(0); circular polarizations for transmitting and receiving to 

reduce transmit-receive coupling, 

where a is the wavelength ratio. Element spacing around the cylinder must be somewhat 

This can be written a.PHI(i)+(dPHI(i)-a.uPHI(i)) and the 25 greater than a half wavelength to avoid the slots from 

term dPHI(i)-a.uPHI(i), which is at least a single constant, running into each other, although it is possible to stagger 

has to be determined in some way to translate the A- or alternate slots by a small vertical displacement to reduce 

C-matrix coefficients determined from receiving mobile their potential mechanical interference or electrical coupling 

signals to the coefficients that shall be used for transmitting with each other. If, for example, 0.75 wavelength spacing is 

to the mobiles. This can, for example, be done by a fixed 30 used, then cylinder circumference is 6 wavelengths, that is 

system calibration that is carried out with the help of a few a cylinder radius of less than one wavelength or about one 

monitoring stations or "dummy mobiles" located at different foot. Such an antenna is considerably smaller than conven- 

positions throughout the service area. Alternatively, by hav- tional three sector antennas. A number of rings of such slots 

ing the mobiles also measure a limited number of residual are stacked vertically with between, for example, 0.5 and 

correlations with signals other than their own, and report 35 one wavelength vertical spacing to provide the same vertical 

these correlations on the slow associated control channel aperture and, therefore vertical directivity, as conventional 

(SACCH), the system can receive enough information to cellular base station antennas. Slots that lie in a vertical 

perform the necessary calibrations for phase mismatch con- column can be connected by feedlines 1801 that feed them 

tinuously. Such reported information can also facilitate in phase. The eight feedlines corresponding to the eight 

calibrating out amplitude mismatches if required. 40 columns of slots are then connected to eight RF processing 

The present invention can also be employed to improve channels 1802. Each RF processing channel comprises a 

the capacity of landbased cellular radiotelephone systems. transmit-receive duplexing filter 1803, a linear transmit 

Such systems generally employ 3-sector antennas to illumi- power amplifier 1804, an RF amplifier 1805, a 

nate three adjacent cells from the same site, as described downconvertor, IF filter, amplifier and A/D converter 1806 

above. Because isolation between sectors is not high (in fact 45 for each frequency channel, and a corresponding transmit 

isolation is almost zero for a mobile on the border of two modulator 1807 for each frequency channel, the outputs of 

sectors), it is not possible with conventional systems to which are summed in summer 1808 before being amplified 

permit use of the same frequency channel in all three sectors. in power amplifiers 1804. 

According to exemplary embodiments of the present The digitized outputs for all eight columns of slots for 
invention, however, the same channel can likely be 50 each frequency channel are fed to a receive matrix processor 
employed as many times as there are antenna elements to 1809. The receive matrix processor 1809 is analogous to the 
form sectors. Thus a three-sector antenna (typically formed matrix processor 650 of FIG. 9. The matrix processor 1809 
by three vertical collinear stacks of dipoles in a corner separates signals arriving on the same frequency but from 
reflector) provides the opportunity to re-use the same chan- different angles such that cochannel interference from 
nel three times. 55 mobiles in communication with the same site is substantially 
Land-based cellular communication capacity is limited by suppressed. The separated signals are fed to voice or random 
the parameter of carrier to co-channel interference ratio access channel processors (not shown in FIG. 18) analogous 
(C/I). The C/I which would be obtained if signals on the to channel processors 660 of FIG. 9. Correlation measure- 
same frequency are radiated around 360 degrees of azimuth ments performed by the channel processors (not shown) are 
is the same as the C/I which would be obtained with 60 fed to a control processor (not shown) analogous to control 
centrally illuminated cells. A 3-cell cluster or site then processor 1702 in FIG. 17. The control processor (not 
becomes the equivalent of a centrally illuminated cell as shown) produces both receive and transmit matrix coeffi- 
regards the re-use pattern needed to achieve a given C/I. It cients for receive matrix processor 1809 and transmit matrix 
is known that a 21 -cell re-use pattern is needed to provide processor 1810 to produce a transmitted signal to every 
the required C/I in the AMPS system, therefore a 21-site 65 cochannel mobile in a non-interfering manner, 
re-use pattern would be needed if all sectors in the same site A difference in propagation conditions can arise in land- 
used the same frequencies over. This compares with the mobile applications as compared to satellite applications, 
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resulting in some modifications to the matrix processing that It is thus more convenient to regard the C-matrices as 

will now be described. Satellite propagation paths are sub- constant relative to directions of arrival, and to introduce an 

stantially line of sight, and even if signal echoes from explicit set of Rayleigb fading variables to explain the fast 

objects in the vicinity of the mobile occur, the sightlines fading. Each signal in the vector Tj, the first signal tl(j) for 

from these objects to the satellite are substantially the same 5 example, thus has an associated complex multiplying factor 

as the direct ray from the mobile to the satellite when f1 q representing the undelayed Rayleigh fading path from 

compared with the relatively large cell diameters in satellite- mobilc i to the array. Assembling the factors rl,r2,r3 . . . m 

cellular systems. doWT1 me diagonal 0 f a ma trix, with zeros elsewhere and 

This is not true for landmobile systems. A substantial echo denotin this fadin matrix b R0 ^ M of ^ si ^ 

from a large building or mountain ranee on the other side of . , . . 

(U & j . 4L i -| u . , lL , 10 are then simply given by: 

the antenna compared to the mobile can result in an echo that r j &- j 

comes in from a direction anywhere between 0 and 180 /ear/ 

degrees away from the direct ray. Since such echoes carry ^ „ . p . ^ ¥ „ , „ 

signal energy, it is often desirable to exploit them to provide Defining a different fading matrix RI for the first delayed 

a diversity path in the event that the direct ray fades or is P ath > me dela y cd faded ^ h are S iveD b ? : 

shadowed in order to improve reception. Typically, the 15 RiT(J-m) 

signal path from the mobile to the base station antenna _ _ t t t 

consists of a number of rays caused by reflections from nus the s, S nals out of mc dements are given by: 

objects close to the mobile; these rays are received substan- Sj<:JtQ. 7J+C' J?i . T(j-m) 

ti ally from the same direction and combine to produce . _ , . 

so-caUed Rayleigh fading. Since the base station antenna in, 20 J f?* a f. 10 °? e T^n °l , ,nvenUo °' ^P^ 00 

f I I 11 i *- • j vu * i i J of the fading signals RO.Ti takes place using the separated 

for example, large-cell applications is deliberately placed . , „ A Jt. v i , i ^ . . t 

high at a good vantage point, there are not expected to be s, 8 na ^ R0T &- m ) calculated m samples ago, based on the 

large reflecting objects in close proximity, for example ct * ua on * 

within 1.5 Km, that could result in rays coming from r , 

substantially different directions. This means that rays 25 RQ Tj = C- 1 - 1 Sj-C- -~- ■ (f(0-T{j- m)) I 

reflected from such objects and coming from an arbitrary * 

direction would be expected to have traversed a larger . . . 

distance, e.g., 3 Km, and thus suffered a delay of 10 /iS or li 15 seen mat me Previously separated signals RO.TO-m) 

must first have their fading factors removed by division by 



more. 



To take care of both types of the aforementioned 30 R0 to re P lace me fadm S factors for the direct ra y s ^ the 

phenomena, that is a cluster of rays from substantially the fadin g faclors R1 for the M*yvd rays. This can cause 

same direction causing the signal to exhibit Rayleigh fading numerical difficulties when a signal fades out completely so 

as well as a cluster of rays from a substantially uncorrelated mal ils associated r-factor becomes zero. However, since the 

direction representing a delayed signal, another term can be separated signal would also become zero, it is possible to 

introduced into the receive matrix processing as follows. 35 a meaningful value to R0.T(j-m)/R0, using, for 

A signal sample Si(t) received at the ith antenna element example, knowledge of the nature of the transmitted signal, 

(column of slots) is the sum of non-relatively-delayed trans- For example, knowledge that the transmitted signal is a 

mitted signals Tk(t) from mobiles k and signals relatively ' onstant amplitude signal, or that it should be continuous 

delayed by dt given by: between samples, could be used. 

40 Alternate implementations and modifications of the prin- 

Sf(i)=Cfl.n(f>+07J2(r) . . . -fCw.7h(f>+c/i'.7i(f-rf/>+c/2+.72(f- ciples set forth herein will be readily appreciated by those 

dt) . . . +Cin'.Tn(t-dt) skilled in the art. For example, although the signal trans- 

When the equations for all Si(t) are collected into matrix mitted in ^ foturc cellular system will most likely be a 

form, they can be written: digital signal, the principles of the present invention are also 

45 applicable to analog signals. In both cases, the fading spectra 
Sj^.Tj+C .TQ-m) ^ e ^ ^ p our j er transform of a successive series of r-values) 
Wherein the suffix j of T means values at a current time and are narrowband compared to the modulation, which is the 
the suffix j-m means values m samples ago, corresponding means by which information in the modulation can be 
to the delay dt. For example, if the signals are sampled every distinguished from modulation caused by fading. In the case 
5 //S, then for a delay dt-10 j/S, m would be equal to 2. 50 of digital signals, the modulators used at the transmitters are 
The signal fading of the undelayed ray can be considered well characterized a-priori, so that the waveforms Tj that 
to be due to varying C coefficients, with the transmitted they will produce for a given information bit pattern can be 
signals T being constant, or the signals T can be considered predicted. If a known bit pattern is contained in a segment 
to be varying due to Rayleigh fading and the matrix C to be of transmission, a corresponding segment of the Tj wave- 
constant. The latter is considered here, because after sepa- 55 forms can be predicted and correlating this with the received 
rating the fading signals T by using constant matrices, the signals will yield an estimate of the corresponding T-value. 
voice channel processors can handle the fading signals as This process is referred to as "channel estimation". The 
they do in landmobile systems. channel estimates may be updated after decoding each 
If the signals T are considered to be fading, however, note information bit. Due to the channel varying much slower 
that the fading on the delayed term is not correlated. In order 60 than information bits and even more slowly than the sample 
to be able to consider T(j-m) as a delayed replica of the rate of Tj, which may be, for example, eight times the 
fading signals Tj therefore, the difference in fading must be information symbol rate, channel estimates are averaged 
explained by regarding the coefficients C as varying to over many successive samples of the T- wave forms, and are 
convert the fading on the direct ray to the fading on the thus somewhat less noisy than the information signals 
delayed ray. However, infinite values of C would then arise 65 themselves, 

due to the varied coefficients being the ratios of Rayleigh In the case of analog FM signals, for example, the 

fading values. modulation is known a-priori to be constant amplitude, 
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varying only in phase. The rate of change of phase is known 
a-priori to be restricted to a value corresponding to the 
maximum frequency deviation, and the frequency variation 
is continuous and so the phase and at least its first and second 
derivatives are continuous. This a-priori knowledge can be 
used to predict a next Tj value from the previous history. For 
example, if Q^- was the previous phase estimate and Q its 
derivative estimate, and was the previous amplitude 
estimate, then Ty-A,. EXP and T^-A, EXP 

(,{Q,+Qdt)). Hence, T >+1 is predicted from T J+1 -T EXP 

Channel estimation techniques often use a Kalman filter 
including derivatives, in which a prediction of the next value 
of the channel estimate is made using an estimate of the time 
rate of change (derivative) of the signal, then the predicted 
channel estimate is used to predict the next signal sample 
point. The error between the predicted and received signal is 
then used to correct the estimate of the channel (the fading 
factor) and its derivative in such a way as to sequentially 
minimize the sum square error. 

The same Kalman filter technique can also be used to 
estimate the diagonal elements of both RO and Rl. Having 
estimated these diagonal values, according to another aspect 
of the present invention, it is ascertained whether any value 
of Rl is greater than a corresponding value of RO. If a value 
of Rl is greater than a corresponding value of RO, that would 
indicate that the delayed ray is currently received at a greater 
strength than the direct ray. Then the column of C corre- 
sponding to that element of Rl is swapped with the corre- 
sponding column of C corresponding to RO to form new 
matrices which are denoted by Cmax and Cmin. The greater 
element from Rl is swapped with the corresponding smaller 
element from RO to form new R-matrices Rmax and Rmin, 
respectively. The elements of T(j-m) corresponding to the 
swapped R-elements are then swapped with the correspond- 
ing elements of Tj to form mixed vectors of delayed and 
undelayed signals denoted by Uj and Vj, respectively. The 
vector Uj can contain some elements of Tj and some 
elements of Tj-m), while the vector Vj then contains the 
remainder. Thus, the equation for signals out of the array 
elements becomes: 

Sj^max Jl max. Uj+CminJZmui.Vj 

This equation can then be solved to yield: 

J//-{CmaxJ?max]" 1 .[5/-CminJ?min.V7] 

Since each element of Rmax was chosen to be the greater of 
two, the chances of zero values are reduced. Furthermore, 
the Vj values that have to be subtracted from Sj are mini- 
mized by multiplication by Cmin, so if Vj values are wrong 
or noisy the error propagation into subsequent values win be 
attenuated. 

The vector Vj, however, contains some as-yet uncal di- 
lated values. Assuming that the same elements of RO and Rl 
are chosen for Rmax and Rmin next time, the as-yet uncal- 
culated values of Vj belong to a future U-vector U(j+m). The 
previously calculated values of T contained in Vj come from 
a previous U-vector, U(j-m). 

Cmin and Rmin can be partitioned into two matrices 
Cminl, Rminl and Cmin2,Rmin2, the columns of which are 
associated with the Vj values that come from previous or 
U-vectors, respectively. Thus, the U-vectors can be 
described as: 

l//=|CinaxJimax]" 1 l5)'-C[ninl Jfminl. £/0'-'")-Cmin2J&nin2. U{j+ 

The values of U(j-m) are known from a previous 
calculation, but the values of U(j+m) are not. Therefore, Uj 
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is first calculated on the assumption that all U(j+m) are zero. 
Then, m samples later when U(j+m) has been calculated on 
the assumption that U(j+2m) are zero, the calculated values 
of U(j+m) can be back-substituted into the above equation to 
give a refined set of values for Uj. These Uj values may be 
then back-substituted into a previous calculation of U(j-m) 
to refine that calculation, and/or forward substituted into the 
calculation of U(j+m), or both, to an iterative extent limited 
only by available processing power in the receive matrix 
processor. 

Simplifying the above equation by denoting: 

Ao-[CinaxT l 

A 1-tCmaxT 1 {Cminl Jim in2] 
A2-[Cmaxy l {Cmin2J2mic2] 

and substituting yields: 

nmax.Uj^o.S]-Al.UO-m)-AZU(j+m) 

If Al has diagonal elements Dl and A2 has diagonal 
elements D2, then we can also write: 

Dl. Ufj-m)+Rm2x. E//+D2. U(j+m)=Ao.Sj-(A l-Dl).U(j-my 
(A2-D2).U(j+m) 

The left hand side of the foregoing equation represents 
separated signals without cancellation of delayed or 
advanced rays. The separate channel processors can process 
these signals including delayed echoes to obtain better 
quality demodulation and decoding than if echoes had been 
subtracted. The improved decoded signals are useful in 
better producing the required channel estimates. A device 
that can, for example, be used for this purpose is a Viterbi 
equalizer such as described in commonly assigned U.S. Pat. 
No. 5,335,250 entitled "Bidirectional Demodulation Method 
and Apparatus", which is hereby incorporated by reference. 

Thus, according to this exemplary embodiment of the 
invention, echoes of each signal are subtracted from esti- 
mates of other signals, but not from the estimate of the signal 
itself, to produce separation of signal+echo signals that are 
processed by individual channel processors. Echoes of each 
signal itself are left in additive combination with the signal 
and are used by a Viterbi equalizer. If echoes are not delayed 
or advanced by multiples of the modulation symbol period, 
a so-called fractional-spaced Viterbi equalizer can be used. 

Such equalizers continuously estimate and update the 
amount and phase of additive echoes, as described in com- 
monly assigned U.S. Pat. No. 5,164,961 to Bjora Gudmund- 
son entitled "A Method and Apparatus for Adapting a Viterbi 
Algorithm to a Channel Having Varying Transmission 
Properties", U.S. Pat. No. 5,204,878 to L. Larsson entitled 
"Method of Effecting Channel Estimation For a Fading 
Channel When Transmitting Symbol Sequences", and U.S. 
Pat. No. 5,465,276 and entitled "A Method of Forming a 
Channel Estimate for a Time Varying Radio Channel", each 
of which are incorporated here by reference. The estimated 
values correspond to the diagonal elements of the diagonal 
matrices Dl, Rmax, D2. Knowing Cmax and Cmin, Rminl 
and Rmin2 can then be determined, thus the channel adap- 
tive equalizers in the individual channel processors can 
determine the Rayleigh fading functions R0 and Rl. 

A purpose for cancelling by subtraction cross-echoes, i.e., 
echoes of one signal that are additive to a different signal, is 
to provide separate signal sample streams that each depend 
only on one signal and its own echoes, as such can be 
handled by said channel-adaptive, Viterbi equalizers. For 
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completeness however, a further method will now be binary bit of each signal. How it makes this postulate is 

explained, that can be used when the number of signals to be immaterial, as all postulates will eventually be tried. In the 

separated is relatively few, for example, eight signals. event that postulates are tried sequentially, they can, for 

The receive matrix processor can be regarded as undoing example, be generated by an 8-bit counter. If however all 

the additive signal mixing that takes place in the aether. This 5 postulates are tried in parallel using replicated hardware, 

is advantageous in simplifying the operation of the channel each hardware unit would handle one fixed postulate that 

processors. However, as disclosed above, numerical diffi- could then just be hard-wired in. 

culties can arise in dealing with signals that can periodically Together which each of the previous 8-bit postulates in 

fade completely. This can result in certain matrices becom- turn plus the new 8-bit postulate, a set of eight signal 

ing singular, i.e., difficult to invert accurately. An equivalent 10 predictors 1930 predicts the complex value of each signal 

problem arises in equalizers that attempt to undo the effect incident on the array including one or more reflected rays by 

of a corruptive channel, for example a channel that suffers using the fading channel coefficients R and R' and a-priori 

from selective fading that causes a null in the transmission knowledge of the transmitted modulation or coding. The 

function at some frequency. An inverse channel filter that complex signal values are then combined in matrix proces- 

attempts to undo the effect of such a channel would try to 15 sor 1940 by calculating the equation: 
create infinite amplification at the null frequency, with Q-^Rr+TR'Tr- 
consequent huge amplification of noise and other difficul- J • J • 0- m ; 

ties. where C and C ! are square matrices representing the direc- 

Therefore it is often proposed, as in the Viterbi equalizers tions from which the direct and delayed waves are princi- 

cited, that the channel should not be "undone" by subjecting 20 pally received. 

the received signal to an inverse channel filter to produce an The calculated signals Sj are the signals that are expected 

undistorted signal that is then compared to the alphabet of to be received at the array elements if the hypothetical eight 

expected symbols, but rather the alphabet of expected sym- bits are correct. These hypothetical signals are then com- 

bols is subjected to the same channel distortion as the signal pared with the corresponding received signals R1,R2 . . . R8 

by use of a mathematical model of the channel, and the 25 from the array elements using comparator 1950. Comparator 

distorted received signal compared with this predistorted 1950 evaluates the net mismatch of the eight predictions 

alphabet. from the eight array signals by, for example, computing a 

According to a further exemplary embodiment of the sum of squares of the differences. Other means to produce 

present invention, a method is disclosed whereby no attempt a signal representative of the net mismatch are however 

is made to separate or "unmix" in the receive matrix 30 known to the art, based on a mathematical expansion of the 

processor the plurality of co-channel signals received by the sum of squares, and can be used if considered advantageous 

array to produce separated signals that are then compared in for the particular implementation chosen. For example, note 

separate channel processors with the alphabets of expected G. Ungerboeck, "Adaptive Maximum Likelihood Receiver 

symbols. Instead, the alphabets of expected symbols are For Carrier Modulated Data Transmission Systems", IEEE 

premixed in every possible way with the aid of a model of 35 Trans. Commun. \bl. COM-22 No. 4, pp. 624-636, May 

the mixing process that takes place in the aether, (i.e., with 1974, U.S. Pat. No. 5,031,193 to Atkinson et al., and U.S. 

the aid of the C-matrix coefficients and the channel estimates Pat. No. 5,191,548 to Backstrom et al., each of which is 

R) and the mixed alphabet is then compared to the mixed incorporated here by reference. The sum square error signal 

signals received by the array elements. is fed back to the SMLSE controller 1910 which adds the 

Such a scheme expands the number of possible mixed 40 error to the previous error stored in state memory 1900 

symbols in the alphabet exponentially according to a power against the previous 8-bit signal hypothesis 1921 employed 

of the number of signals. For example, suppose each signal in signal predictors 1930 to produce the signals ri'ti'. 
is modulated with binary symbols. The expected symbol The above procedure is carried out for each new 8-bit 

alphabet has only two symbols, 0 or 1. However, if the array hypothesis in turn preceded by each of the stored, previous 

elements receive weighted sums of eight signals, each of 45 hypotheses. This results, for each new hypothesis, in 256 

which instantaneously may be modulated with a 1 or a 0, the candidate cumulative error numbers depending on which 

number of possible mixed signals that can be received is 2 8 preceding hypothesis was used. The lowest of these is 

or 256, if all the signals are aligned in time. If different selected to become the new cumulative error associated with 

signals are not time aligned, then a symbol period of one the state corresponding to the new 8-bit hypothesis. When 

signal may overlap two symbols of another signal. Thus the 50 all possibilities for the new 8-bit hypothesis have been 

waveform over a symbol period of one signal can depend on processed in this way, the state memory 1900 will contain 

two symbols of each of the other signals. Nevertheless, each 256 new cumulative error numbers associated with each new 

point of the waveform depends only on the one symbol of hypothesis, as well as a record of the best preceding hypoth- 

each signal whose symbol period it lies in. When echoes are esis to each, i.e., that giving the lowest error, and the 

taken into account, however, each waveform point can be 55 preceding hypotheses to those in turn, and so on. Thus each 

dependent on two symbols of each signal thereby raising the of the 256 states contains a candidate demodulated sequence 

number of possible values that can be observed to 65536. It of 8-bit values. The oldest values in these sequences will 

will however be described below how, for example, a tend to agree and when this happens the machine is said to 

256-state Viterbi algorithm can be used to jointly demodu- have converged to an unambiguous decision. The decided 

late the signals from the array. 60 8-bits are then extracted to yield one bit decision for each of 

According to an exemplary embodiment aspect of the the eight incident signals. If convergence does not occur and 

invention, and referring to FIG. 19, a numerical machine has the sequence memory 1900 becomes full, the path history is 

256 sets of memory banks 1900 each associated with a truncated by believing the oldest byte of the state having the 

specific 8-bit postulate for one previous binary bit in each of lowest cumulative error. That value is then extracted and the 

the eight signals, on which, due to a delayed echo, the 65 path history memories shortened by one byte, 
received array signals will depend. The SMLSE controller The above process represents an alternative to attempting 

1910 now makes another 8-bit postulate 1920 for the current to separate signals that have been mixed by means of matrix 
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processing. Instead, signals are hypothesized by models of 
the transmitters and models of the mixing process, and the 
hypothesis best corresponding to the observed, mixed sig- 
nals is determined by the SMLSE machine 1910 in the 
manner described above. Thus, the need to invert a mixing 
process to separate mixed signals, which may be mathemati- 
cally intractable, is avoided by instead applying the math- 
ematically tractable mixing process to the hypothesized 
signals to predict the mixed signals that should be received 
by the array elements and picking the hypothesis that best 
matches the observed signals. This process will not fail 
when two mobiles using the same channel lie at the same 
bearing, the process then being equivalent to joint- 
demodulation as, for example, disclosed in U.S. Pat. No. 
5,506,861 entitled "A Method and Apparatus for Joint 
Demodulation of CDMA signals with Multipath Time Dis- 
persion". 

The above-described exemplary embodiments of the 
present invention are applicable to satellite cellular commu- 
nications systems to provide greater use of available band- 
width by permitting immediate spectrum re-use in adjacent 
cells. These techniques have also been described in relation 
to land cellular systems, where they permit, for example, 
re-use of the same frequency in adjacent sectors. 

Ia practice, Inn both the satellite and land-based applica- 
tions of the present invention, benefits are achieved by a 
combination of adaptive signal processing techniques linked 
to traffic management techniques. The traffic management 
techniques relate to continuously operational systems using 
TDMAor FDMAor a combination thereof in which calls are 
continually being terminated and new calls established. By 
selectively establishing new calls on time- or frequency- 
slots in such a way as to optimize a communications 
criterion, a natural sorting of traffic into groups using the 
same timeslot and/or frequency is established. The criterion 
relates to the ease with which the adaptive signal processing 
can separate signals on the same frequency and/or timeslot 
based on the reception of different, linearly independent 
combinations of them using a plurality of antenna elements. 

According to yet another exemplary embodiment of the 
present invention, the signal processing does not adapt to the 
movement of mobile phones or to new call set-up and 
termination, but operates in a deterministic way and instead 
the traffic is adapted to the deterministic characteristics of 
the signal processing using a dynamic traffic channel assign- 
ment algorithm. 

Conventional land-based cellular systems typically 
employ so-called sectorization, in which a single antenna 
mast carries three, 120-degree coverage antennas and illu- 
minates three cells from a common site. This saves on real 
estate costs compared to illuminating the three cells using 
three separate antenna sites at the cell centers. Six sector 
systems are also known. Cellular systems have convention- 
ally employed analog FM voice transmission in which each 
conversation is assigned a separate pair of up-and down-link 
frequency channels respectively. TO MA systems are now 
being installed using digital speech transmission, in which 
each conversation is allocated a unique pair of timeslot- 
frequency channel combinations. In these conventional 
systems, however, the three, 120-degree sector antennas 
have the same radiation patterns for all frequencies and/or 
timeslots. 

According to yet another exemplary embodiment of the 
present invention, rotationally offset radiation is provided 
between different frequencies and/or timeslots. For example, 
on frequency channel 1 the three, 120 degree sectors may be 
orientated towards 0 degrees (Due North), +120 degrees 
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(South East) and +240 degrees (South West). On frequency 
channel 2, the three sectors may be orientated to 60 degrees 
(North East), 180 degrees (Due South) and 300 degrees 
(North West). In general one might have as many as 120 

5 frequency channels with corresponding antenna sector pat- 
terns offset by only one degree from each other. Such a 
system cannot be implemented using today's fixed-beam 
cellular base station antennas, but can be arranged with 
using the exemplary cylindrically symmetric array and asso- 

10 ciated matrix processing of FIG. 18. 

Similarly, the antenna sector patterns can be rotationally 
staggered as between different timeslots in a TD MA system. 
In either the FDMA or TDMA or hybrid cases, this exem- 
plary system determines at call set up, and optionally at 

15 regular intervals thereafter, the optimum time- and or fre- 
quency slot combination to use for communicating with the 
mobile station. The combination of a frequency and timeslot 
is abbreviated henceforth to simply "channel". The optimum 
channel is most likely one which has an associated antenna 

20 sector pattern pointing in the direction of that mobile. This 
channel would be selected if the selection criterion is, for 
example, maximum signal strength and the channel was 
free. If the criterion is maximum signal-to-interference ratio, 
different selections can result. Adaptive channel selection 

25 methods can be used to implement the present invention as, 
for example, disclosed in U.S. PaL No. 5,230,082 to Ghisler 
et al. which is incorporated here by reference. 

FIG. 20 illustrates a set of staggered sector patterns that 
can be produced by the arrangement of FIG. 18 using fixed 

30 matrix coefficients for each frequency channel of an FDMA 
system. Three lobes are created in this example on every 
frequency channel. The notation Pi(Fk) indicates the pattern 
of the i lobe on the k** frequency channel. The matrix 
processing coefficients are preferably chosen such that 

35 Pl(Fk) and P2(Fk) have minima where P3(Fk) has its 
maximum, and reciprocally. If the minima are zero, the three 
lobes are said to be orthogonal. That permits a mobile 
located in the nulls of PI and P2 to receive maximum signal 
from P3 with no interference from the other two, which can 

40 thus carry separate signals. In general, true zeros will not be 
perfectly achieved, and the channel selection criterion will 
thus allocate a mobile to a frequency where the correspond- 
ing sector patterns result in maximum ratio of wanted signal 
to unwanted interference from other lobes and other cells. 

45 For example, the mobile M in FIG. 20 would be allocated 
preferably F4, where the lobe P3(F4) has the maximum 
strength in the direction of the mobile M. If P3(F4) was not 
available, the next best allocation P3(F3) would be tried, and 
so on. 

50 In practice, an FDMA cellular system such as AMPS has 
1000 channels available, usually divided between two 
operators that handle a minimum of 400 each. Using the 
traditional frequency re-use pattern of 21, this results in 
around 20 frequencies being available in every cell or 60 per 

55 site. The angular difference between lobes on different 
frequencies would thus, in a three lobe system, be only V&oth 
of 120 degrees or 6 degrees. In this example, different lobes 
at the same site all have different frequencies. Assuming 
uniform distribution of mobiles in angle, the channel allo- 

60 cation algorithm would result in each mobile being within a 
few degrees of beam center. This results in mobiles receiv- 
ing better signals on average than in today's fixed sector- 
ization patterns which, when optimized, are around 12 dB 
down at sector edges. If the wanted signal is improved in this 

65 way, the tolerance of interference from surrounding cells is 
improved such that the re -use pattern can be shrunk from 21 
to a tighter re-use pattern such as 12, with a consequent 
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capacity gain of 21/12. This can be achieved using the same 
number of lobes as sectors in today's cellular systems. If the 
number of lobes is increased to eight, as illustrated in FIG. 
18, a further 8/3 increase in capacity is obtained, to around 
five times current AMPS capacity. Moreover, allowing every 
cell to adoptively select any of the 400 frequency channels 
in attempting to maximize signal to interference ratio gains 
a factor of two in capacity relative to having a fixed subset 
of frequencies (V&ist or l /i2th of the total) in each cell. This 
is achieved when transmit power levels are also adapted to 
the varying radial distance of each mobile from its cell site. 
It is also possible to use all 60 site frequencies in each 120 
degree sector by making lobes using the same frequency 
orthogonal, as defined above. Lobe separation is then 2 
degrees and the channel allocation algorithm ensures not 
only that each mobile is within a couple of degrees of beam 
center, but also within a couple of degrees of the minima of 
the co-frequency lobes. 

If instead of associating staggered sector radiation pat- 
terns with different frequency channels Fl JF2 ,F3 . . . they are 
associated with different timeslots of a TDMA signal using 
a single frequency, the resulting radiation from the base 
station antenna will take a certain set of directions for 
timeslot 1, a set of rotated directions for timeslot 2 and so 
forth, such that the beams are apparently rotating with time. 
Thus in the TDMA context this exemplary embodiment of 
the present invention may be formulated in terms of creating 
beams which continuously rotate through 360 degrees over 
a TDMA frame, or more appropriately, rotate by 360/N 
degrees during a TDMA frame where N is the number of 
sectors of frequency re-use, and the data modulation for the 
next frame is shifted back one sector between successive 
frames such that data for the same mobile continues to be 
radiated in the same direction. Data destined for a particular 
mobile is indicated in the US IS54 TDMA system by 
inclusion of a "Digital Voice Color Code" (DVCQ in 
TDMA bursts. Thus, for example, this technique can be 
described more simply in terms of rotating the antenna 
sector patterns in one direction while rotating the DVCC in 
the reverse direction at the same rate such that the same 
DVCC continues to be radiated in the same direction on 
successive frames. 

Both exemplary FDMA and TDMA embodiments of the 
present invention provide mobile stations with the capability 
to determine coarse geographic position. In the FDMA 
version, the mobile measures relative signal strength on 
different frequencies. The frequency on which the greatest 
signal strength is received indicates the bearing of the 
mobile within a sector. The sector is determined by decoding 
sector ID information contained within the transmission. 

In exemplary TDMA embodiments, the mobile does not 
even have to change frequency. The mobile instead notes the 
cyclic signal strength variation during a TDMA frame and 
then determines the peak and trough signal strength posi- 
tions relative to timeslot 1, which can be identified by the 
slot ID information carried in each slot. 

The cyclic signal strength variation can be processed over 
several cycles with the aid of a Fourier transform and the 
phase of the fundamental component relative to timeslot 1 
will then indicate the mobile's bearing. Bearings from two 
base stations of known positions then fix the mobile posi- 
tion. The mobile can report the timing of signal strength 
peaks and the network can perform the position calculation, 
rather than the network having to send coordinates of base 
stations to the mobiles. Upon allocating a trafiBc channel to 
a called or calling mobile, the network can then determine 
the best of all available timeslot/frequency combinations to 
use. 
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The above-described can also be adapted to provide 
advantageous communications between mobile stations and 
an orbiting satellite. According to this embodiment, the 
antenna array signal processing is not adapted to various 

5 mobile positions, but rather mobiles are allocated to a 
specific antenna array signal processing channel based on 
position in such a way as to optimize communications. That 
is, mobiles are adaptively allocated to communicate using 
one of a number of fixed, staggered antenna beam positions 

10 instead of adaptively steering the antenna beams onto the 
mobile positions. 

The operation for satellite use may be modified slightly. 
The notion of fixed antenna beams would be applicable to a 
geostationary satellite, but may not be applicable to, for 

15 example, a low-orbit satellite that changes position relative 
to the earth. Then the position of a beam relative to a given 
mobile would move due to satellite motion if not due to 
mobile movement. If the satellite beams move over the earth 
relatively slowly in comparison with the average 3-minute 

20 call duration, it may be sufficient to allocate a mobile to a 
beam at call set up, as one would in the Geostationary case. 
However, according to this examplary embodiment of the 
present invention, beam directions can be adapted to remove 
the systematic motion of the satellite over the earth so that 

25 the area illuminated by each beam is static from satellite rise 
to satellite set. In this way, a mobile may remain allocated 
to the same beam irrespective of satellite movement during 
this period. 

Furthermore, such a system of low-orbiting satellites 

30 would generally be arranged to provide continuous coverage 
whereby as one satellite sets, another rises. For example, it 
can be arranged that a satellite rising in the west takes over 
the illumination of the same area just being vacated by a 
satellite setting in the east. Then, as the adjacent area to the 

35 east of this experiences loss of the setting satellite, the rising 
satellite creates a new beam to dovetail in while having 
maintained the first beam over the original area, and so on 
until the new satellite has taken over ttlumination of all areas 
originally served by the setting satellite. 

40 Thus, the application of this exemplary embodiment of 
the present^ invention in the case of moving satellites allows 
the illumination patterns from the satellite antenna to be 
compensated for satellite motion so as to illuminate fixed 
areas of the earth while adaptively allocating mobiles to 

45 illumination patterns using a channel assignment algorithm 
that optimizes a communications quality criterion. This 
contrasts with mechanical methods of compensating for 
satellite motion by tilting the satellite or antenna so as to 
maintain the center point of at least one area constant This 

50 mechanical method, however, cannot maintain the center 
points of cell illumination areas constant due to these areas 
changing shape from circular as the satellite moves overhead 
to elliptic and finally to parabolic at earth-edge illumination. 
Alternatively, this exemplary embodiment of the present 

55 invention can employ both a mechanical method for coarse 
compensation plus the method of adaptive antenna array 
signal processing to correct the illumination patterns for 
shape change as the satellite moves. Alternatively, signal 
processing can be used to hold the areas served by a 

60 particular frequency and/or timeslot constant while progres- 
sively creating new areas forward of the satellite's ground 
track that are being vacated by a setting satellite and while 
terminating illumination of areas to the rear of its ground 
track that are being taken over by a rising satellite. 

65 The operation of this exemplary embodiment of the 
present invention is depicted in FIGS. 21(a) and 21 (b). At a 
certain time T (FIG. 21(c)) a rising satellite 2100 illuminates 
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areas with frequencies (left to right) F1,F2 J*3J?1,F2,F$JF1, It will be apparent to one skilled in the art that TDMAand 

F2 and a falling satellite 2102 illuminates further areas with FDMA are not the only access methods that are compatible 

frequencies (left to right) F3,F1,F2,F3,F1,F2,F3,F1 which with the present invention. Code Division Multiple Access 

continue the frequency re-use sequence. At, for example, (CDMA) can also be used, where illumination areas are 

time T+5 minutes (FIG. 21(6)) the rising satellite 2100 has 5 similarly staggered over the earth according to a CDMA 

ceased to illuminate the rearmost Fl area 2104 which is code use pattern. Indeed, any multiple access method which 

presumably now obscured from view (i.e., the satellite is too defines a channel by means of a set of access parameters can 

low on the horizon for good communications with this area) have systematically staggered illumination areas depending 

while the setting satellite 2102 has stopped illuminating its on those access parameters. Moreover, the access method 

rearmost area 2106 with frequency F3 for the same reason. 10 used on the downlink can be different from the access 

On the other hand, the rising satellite has created a new method used on the uplink, providing a set of uplink access 

illumination area 2107 forward of its ground track to fill in parameters is paired with a corresponding set of downlink 

the area vacated by the setting satellite. The rising satellite parameters in each offset beam or staggered illuminated 

2100 can appropriately illuminate the new illumination area area. For example, a combination of TD MA on the downlink 

2107 with the same frequency as its predecessor used. 15 with CDMA on the uplink, the uplink transmission being 
Meanwhile, the satellite 2102 that is setting with respect to continuous apart from a short interruption during reception 
this area is rising as viewed from areas ahead of its ground of the downlink slot, is disclosed in U.S. patent application 
track, and uses the released capacity to create a new area Ser. No. 08/179,954 filed on Jan. 11, 1994 and en tided 

2108 forward of its ground track illuminated with frequency "TDMA/FDMA/CDMA Hybrid Radio Access Methods", 
F2. 20 which is incorporated here by reference. Having described 

It will be appreciated that instead of different frequencies, an exemplary embodiment wherein dynamic traffic channel 

the overlapping areas could have been allocated different assignment allows traffic to be adapted to the deterministic 

timeslots in a TDMA frame, or different frequency/dmeslot characteristics of signal processing, a complementary, 

combinations in a hybrid FDMA/TDMA system. Either way, exemplary embodiment will now be described wherein 

the availability of a large number of channels allows the 25 capacity can be optimized through coding and frequency 

overlapping beams to be much more finely spaced than in re-use schemes. 

the example of FIGS. 21(a) and 21(b), so that it is almost The ultimate capacity of a cellular satellite communica- 
equally effective to allocate a mobile to an adjacent beam as tions is limited by available bandwidth, as power limitations 
to the optimum beam. Logically, one should allocate a can always be solved by money, e.g., by launching more 
mobile preferably to a be am in which the mobile is centrally 30 satellites. Practically, however, there are financial con- 
located. However if the corresponding channel is occupied, straints on power and political constraints on bandwidth, 
the mobile can be allocated to a slightly off-center beam and therefore it is desirable to use bandwidth efficiently without 
may be handed over to the on-center beam when the call significant sacrifice of power efficiency, 
using that channel terminates. It shall be appreciated that the trade-off of bandwidth and 

In the exemplary TDMA embodiment, the rising satellite 35 power efficiency for a cellular (i.e., area or global coverage 

and the setting satellite can both illuminate the same area system) is different from that of a single link, as a single link 

using the same frequency, providing different timeslots are trade-off does not consider the possibility of frequency 

used. Thus, a channel and satellite allocation strategy re-use in adjacent cells. The units of capacity in the two 

according to the present invention is to allow calls in the cases are in fact different, being Erlangs/MHz for a single 

changeover region to terminate naturally on the setting 40 link and Erlangs/MHz/SqKm for a cellular system, 

satellite and to re-employ their vacated timeslots in the same A cellular system illuminates a service area by dividing it 

region and on the same frequency to set up new calls using into cells and using some fraction 1/N of the total available 

the rising satellite. bandwidth in each. A cluster of N neighboring cells can thus 

FIG. 22 is a block diagram of an exemplary control be allocated different 1/N fractions so that they do not 

processor that supplies the matrix coefficients to the numeri- 45 interfere. Outside the cluster, for cells far enough away, the 

cal matrix processor of the hub station, e.g., block 1603 in bandwidth can be re-allocated to another cluster. 

FIG. 13. Inputs to the control processor 2200 include The reduction of interference by employing an N-cell 

satellite orbital data including attitude control information re-use pattern is measured in terms of the carrier to inter- 

from the independent satellite Telemetry, Tracking and ference ratio C/I, which is the ratio of wanted signal power 

Command (TT&Q subsystem (not shown). Using satellite 50 to the sum of the power of all unwanted spectrally and 

orbital and satellite antenna pointing information (attitude temporally overlapping signals. Increasing N increases the 

control information) and an input from a real time clock, the C/I, but reduces the bandwidth available in every cell, 

control processor 2200 can determine the matrix coefficients thereby limiting system capacity. Reducing N worsens C/I 

needed such that a given area will be illuminated by a but increases the bandwidth available to every cell. If the 

specific frequency in a specific time slot These coefficients 55 modulation and coding scheme can tolerate the reduced C/I, 

are systematically updated in step with changes in the real capacity will thus be increased by reducing N. 

time clock to maintain these illuminated areas approxi- One method of providing greater C/I tolerance is to use 

mately fixed irrespective of satellite motion. The control redundant coding. This method increases the bandwidth per 

processor 2200 also receives information transmitted from signal, however, which offsets the benefit conferred by 

mobile stations making a random access on a calling channel 60 shrinking the re-use pattern N. The question to be asked is 

that allows the control processor to determine the best where the optimum hes. 

available channel/beam combination to use. This informa- In land-based cellular systems, this question has been 
tion provides a rough indication of the location of the mobile deeply studied, leading some people to conclude that the 
and the control processor then determines the available extreme bandwidth expansion of CDMA techniques corn- 
beam centered most closely on this location. This in turn 65 bined with immediate frequency re -use in every cell pro- 
determines the frequency and/or timeslot that should be used vides the highest capacity. According to exemplary embodi- 
for communication with the mobile. ments of the present invention, however, it is found that 
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capacity increases with increasing coding and reduction of N 
until N=l is reached with a coding rate of about x h (far 
landcellular). At this point the system is not regarded as 
being truly CDMA as each channel is still only used once in 
every cell. CDMA can be denned as the use of each channel 
more than once in each cell, i.e., a fractional value of N. For 
example, N=V£ means each channel is used twice in every 
cell, which would be classified as CDMA. 

Whether this further reduction of N to fractional values 
continues to increase capacity depends on what type of 
CDMA system is employed and on the nature of the propa- 
gation channel and receiver complexity used in the system. 

Three types of CDMA systems may be distinguished: 

i) Conventional, non-orthogonal CDMA 

ii) Orthogonal CDMA 

iii) Interference cancellation CDMA (subtractive CDMA, 
joint demodulation, etc.) 

For the landbased cellular world, it is found that the 
capacity drops off below N=l for CDMA of type (i), levels 
off for orthogonal CDMA (which is really equivalent to 
giving every signal a unique frequency or timeslot) and 
increases for systems of type (iii). Moreover, the gain found 
in systems of type (iii) for landbased cellular where N<1 is 
due to the high near-far environment such that the interfer- 
ence averaging inherent in CDMA techniques includes many 
transmitters of significantly reduced power, and due to the 
landbased cellular scenario being C/I limited rather than 
noise, C/No, limited. Neither of these features are relevant 
to satellite communications systems. Accordingly, the 
present invention explores what kind of coding/frequency 
re-use trade-off would maximize capacity for a given band- 
width allocation in C/No-limited satellite communication 
systems. 

The signal spillover between cells in a landbased cellular 
system is a function of the fourth-power-of-distance propa- 
gation law. In cellular-satellite systems C/I is a function of 
antenna beam pattern sidelobes. It is necessary therefore to 
develop some model of antenna beam patterns to perform 
coding optimization. 

The beam pattern of the antenna depends on the surface 
current distribution over the aperture, called the aperture 
illumination function. Without invoking the supergain 
phenomenon, the most efficient use of aperture is obtained 
with uniform illumination. This gives the best gain but the 
highest sidelobes. The radiation pattern is plotted in FIG. 23 
for a uniformly illuminated circular aperture. The sidelobes 
in the E and H planes are slightly different owing to an extra 
cosine factor that appears in the plane containing the surface 
current vector. This difference manifests itself as cross- 
polarization components when circular polarization is 
employed. Henceforth the E and H plane patterns will 
simply be averaged for the calculation of C/I. 

Reference is made again to FIG. 5 which illustrates a 
3-cell frequency re-use pattern, wherein the shaded cells use 
the same channel fl while the others use £2 or f3. This re-use 
pattern will be used to investigate the coding/frequency 
re-use tradeoff for satellite communications, however, those 
skilled in the art will appreciate that any re-use pattern, e.g., 
7, 9, 12, 21, etc., could be used. Interfering cells lie on the 
points of a hexagon and it suffices to consider the first two 
rings of six interferers. Before their interference levels can 
be computed however, it is necessary to choose the correct 
scaling of the beam pattern to match the cell diameter. 

If the beams are scaled to cross at -3 dB relative to peak 
gain midway between two cells, it is well known that this 
does not result in maximum beam-edge gain. A higher gain 
is achieved if the beam is narrowed, which increases the 
peak gain more than increased edge-loss experienced. 
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FIG. 24 is a plot of peak gain 2401 (at the center of a cell), 
edge gain 2402 (midway between two cells) and the gain 
midway between three cells 2403 as a function of the 
two-cell crossover point in dB down from peak, relative to 

5 the peak gain of the -3 dB crossover case. For reasons that 
will be explained later, the two-cell edge gain has been 
scaled by a factor of two in this plot (i.e., 3.01 dB is added) 
and the 3-cell edge gain has been scaled by a factor of 3 (i.e., 
4.771 dB has been added). This does not affect where the 

lQ respective gains peak, but affects the perception of which of 
the three is the worst case. According to this graph, the worst 
case occurs midway between two cells, and the worst case 
gain is maximized when the 2-cell edge is 3.8 dB down on 
the peak gain, i.e., at point 2404. 

The way in which the C/I parameter depends on the beam 

15 crossover point is shown for the 3-cell re-use pattern of FIG. 
5 in FIG. 25. FIG. 25 is plotted as a function of mobile 
station distance from beam center for crossover points of -3, 
-3.5 and -4 dB showing that C/I over most of the cell radius 
increases if the beams are narrowed beyond that which gives 

20 maximum edge gain. If necessary, choosing a crossover 
point of -4.5 dB would cause negligible loss of edge gain 
while improving C/I at cell center by a further 3 dB to about 
20 dB. C/I at cell edge according to FIG. 4 would be just less 
than 10 dB, but this includes the unmitigated beam edge 

25 crossover loss which, as will be explained later, will not be 
incurred because no mobiles need be located there. 

If mobiles assigned to a particular channel and beam are 
chosen to be those located within 25% of the maximum cell 
radius, the C/l for all points within that area will be as 

30 plotted in FIG. 26. The worst case C/I is maximized to about 
23 dB with a beam crossover design point of -5.5 dB, 
somewhat beyond that which gives maximum edge gain, so 
in practice the -4.5 dB crossover point would be used, 
giving a worst case C/I of 18 dB. 

35 The same calculations are now repeated for the N=l 
frequency re-use pattern, i.e., immediate frequency re-use in 
adjacent cells, and results are plotted in FIG. 27. This shows 
a cell-center C/I of 14 dB for the -4 dB crossover case, but 
a cell edge C/I of about -1.5 dB. The thickness of the curves 

40 in FIG. 27 is due to superimposition of plots for all mobile 
angular positions in the cell, and dependence on angular 
position is a little more noticeable in the N=l case than the 
N=3 case. The first 6 rings of 6 interferers were summed to 
obtain the plot of FIG. 27. 

45 Again, as will be shown later, mobiles using a particular 
channel and beam can be restricted to those lying within 
25% of beam center or less, so it is of interest to maximize 
the worst case C/I within this region, as shown in FIG. 28. 
The worst case C/I in this case maximizes at 13 dB by 

50 choosing the beam crossover point to be -4.8 dB, but this 
may be restricted to -4.5 dB to avoid loss of beam edge gain 
for only a small reduction in C/I to -12.5 dB. 

FIGS. 29-34 give the results for repeating the whole 
process described above for a different aperture illumination 

55 function, the V^-cosine wave. This aperture illumination 
function is slightly less aperture-efficient than a uniform 
distribution, but gives lower sidelobe levels (see FIG. 29) 
leading to higher C/I, particularly in the 3-cell re-use case 
(20 dB over whole cell, or 27 dB out to 25% of radius). As 

60 seen in FIG. 34, the C/I for immediate frequency re-use out 
of 25% of cell radius is 13.5 dB with the practical beam-edge 
crossover point of -4.5 dB. Since this was 12.5 dB for the 
uniform illumination, case of FIG. 28 it should be noted that 
this value is not very sensitive to the aperture function being 

65 used. 

The bit error rate is generally plotted as a function of 
Eb/No, which is equal to the ratio of signal power to the 
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noise power if it were to be measured in a bandwidth equal 
to the bit rate. The latter does not imply any assumption that 
any physical receiver filter bandwidth must be equal to the 
bitrate; it is only that "bitrate" is a convenient unit of 
bandwidth for defining the noise density with which the 
performance of any given receiver will be tested. The 
receiver error rate performance will of course depend on the 
choice of its bandwidth, and that which optimizes perfor- 
mance at a given Eb/No may be greater or less than the 
bitrate depending on the modulation and coding being used. 

The C/I parameter is, by contrast, the ratio of wanted to 
unwanted signal power in the physical receiver bandwidth. 
This ratio, however, is independent of the choice of receive 
filter if the C and I have the same spectral shape and are thus 
equally affected by the filter. With the simplification that any 
T passed through the receive filter will have the same effect 
on error rate as an equivalent amount of white noise, NoB, 
passed by the filter, where B is the noise bandwidth, the 
effect of I can be expressed in terms of an equivalent 
increase in noise density No by an amount Io to No+Io, 
where Io is given by 

I = IoB i.e. Io = IIB 

For BPSK modulation, the optimum receiver bandwidth 
is indeed equal to the bitrate, while for QPSK modulation 
the optimum receiver bandwidth is equal to half the bitrate. 
The bitrate here though is the coded bitrate/chiprate, 
whereas the bitrate for defining Eb/No is the information 
rate. Thus: 

B-Bitrate/r for coding rate r in the BPSK case, 
B=Bitratey2r in the QPSK case, and for general M-ary 
modulation, 

B=Bitrate/rLog 2 (M)=Bitrate/mr where m is the bits per 
symbol. 

Therefore the total bit energy Eb to noise plus interference 
density ratio is given by: 



Eb 

(No+Io) 



[ No , Bitr ate 



For less than 0.5 dB degradation of the Eb/No required for 
a given error rate due to finite C/I, the value of mr.I/C should 
thus be one tenth No/Eb. 

For example, if the ratio Eb/No without interference of 3 
dB is desired, then to operate at 3.5 dB Eb/No the C/I must 
be 10 mr.Eb/No. For BPSK or QPSK and different code 
rates, the required C/I for exemplary coding rates is given 
below: 



REQUIRED C/I using 


BPSK 


QPSK 


Coding rate 1 (none) 


13.5 dB 


16.5 dB 


1/2 


10.5 dB 


13.5 dB 


1/3 


8.7 dB 


11.7 dB 


1/4 


7.5 dB 


10.5 dB 



10 



25 



30 



35 



40 



45 
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The above is for a static channel and does not take into 
account that lower Eb/Nos are needed with lower rate codes 
for the same error rate. 

"Error Correction Coding for Digital Communications" 
by Clark and Cain gives the required Eb/Nos for 0.1% BER 
for constraint length 6 convolutional code rates of 1, %, Yi, 
Vi and V* as follows: 
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r 


Eb/No for BER - 0.1% 


1 


6.7 dB 


3/4 


3.9 dB 


2/3 


3.5 dB 


1/2 


3.0 dB 


1/3 


2.6 dB 



By extrapolation it may be estimated that rate Vi would 
require 2.3 dB with diminishing returns thereafter. Using 
these Eb/No figures, the C/I required for less than a given 
degradation are calculated to be: 



15 



REQUIRED C/I 





for 0.5 dB 


1.0 dB loss 


3.0 dB 


loss 




BPSK 


QPSK 


BPSK 


QPSK 


BPSK 


QPSK 


Coding 


17.2 dB 


20.2 dB 


13.7 


16.7 


9.7 


12.7 


rate 1 














(none) 














3/4 


13.2 


16.2 


10.9 


13.9 


6.9 


9.9 


2/3 


12.2 


15.2 


8.7 


11.7 


4.7 


7.7 


1/2 


10.5 


13.5 


7.0 


10.0 


3.0 


6.0 


V3 


8.3 


11.3 


4.8 


7.8 


0.8 


3.8 


1/4 


6.8 


9.8 


3.3 


6.3 


-0.7 


12.3 


1/5 


5.7 


8.7 


2.2 


5.2 


-1.8 


1.2 dB 



Thus, while the Eb/No for a given error rate levels out 
with increasing coding, the C/I required continues to 
decrease due to the continually increasing bandwidth. This 
equates to the separate concepts of coding gain (which 
applies to Eb/No) and processing gain ( which applies to 
C/I). Coding gain is bounded by Shannon's limit, while 
processing gain continues to increase with bandwidth as in 
a CDMA system. 

The above results for the static channel are pessimistic for 
fading channels. When Rician or Rayleigh fading is present, 
the mean Eb/No can be increased above the static Eb/No 
requirement to maintain the same error rate. However, on the 
satellite downlink the C/I does not exhibit fading, because 
both the I and C reach a given mobile over exactly the same 
channel and fade by exactly equal amounts. Thus the C/I 
does not decrease by 10 dB when the Eb/No fades 10 dB, but 
instead stays at the original value. 

In the fading channel, since the error rate at the mean 
Eb/No is considerably less than the target value, and when 
it fades to the static Eb/No value still only equals the target 
value, it is clear that the error rate only reaches the target 
value by virtue of fades to below the static Eb/No value. In 
fact it can be shown that the preponderance of errors arise 
from the instantaneous Eb/No region well below the static 
Eb/No value, where the same C/I causes less degradation. It 
can be assumed that lower C/I values can be tolerated in 
conjunction with the higher Eb/No values needed to account 
for fading. 

Thus the 12.5-13.5 dB C/I values achievable out to 25% 
of beam radius with immediate frequency re-use are accept- 
able with coding rates of Vz to Vs and using QPSK. Increasing 
the re-use pattern to N=3 would yield 3 times less bandwidth 
per cell, requiring that all coding be removed and even 
higher order modulations than QPSK to be contemplated in 
order to achieve the same bandwidth efficiency, but with the 
penalty of needing considerably higher power (e.g., Eb/No 
of 7.7 dB to 10.7 dB for achievable C/Ts with the 3-cell 
re -use pattern). Thus there is no gain in bandwidth efficiency 
using an N=3 or greater frequency re-use pattern instead of 
N=l, only a major penalty that either is paid in power 
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efficiency (for maintaining the same bandwidth efficiency by 
removing coding, as in the AMSC system) or in bandwidth 
efficiency if coding is retained. 

To make use of the above result it is explained below how 
use of a beam can be restricted to mobiles located, for 
example, only out to 25% of the beam radius or less. 

Again of up to 2:1 in capacity can be achieved in cellular 
systems with a type of frequency planning known as "re-use 
partitioning" In a simple form of re -use partitioning, the 
available channels in a cell are partitioned into three sets that 
are preferentially used for a) mobiles within the inner l A of 
the total cell area; b) for mobiles between l A and *A of the 
cell area, and c) for mobiles in the outer l A of the cell area. 
Assuming a uniform area distribution of mobiles within the 
cell, this partitioning achieves equal demand for each of the 
channel sets. The allocation of channels to equal area rings 
is then permuted in the neighboring cells according to a 
3-cell re-use pattern with the result that no two neighboring 
cells use the same channel out to their mutual border, with 
consequent increase in C/I for no loss of capacity. The 
overall re-use pattern to achieve a given C/I may then be 
shrunk to achieve an increase in capacity. Based on the 
foregoing principles, re-use partitioning and coding can be 
optimally combined according to examplary embodiments 
of the present invention which will now be described. 

FIG. 35 shows a simplified example form for the case 
where three channels or groups of channels (which may be 
frequencies, timeslots, codes or combinations thereof), des- 
ignated by the colors black, red and green are available in 
every beam. The beam edges at the design crossover point 
(e.g., 4.5 dB) are shown by the larger colored circles of FIG. 
35. 

The large black touching circles thus refer to beams using 
the "black" channel and touch at -4.5 dB down from the 
peak gain. The large red touching circles represent the beam 
patterns for the red channel. These are displaced relative to 
the "black" beams, and this fixed displacement is achieved 
for example by modifying the phases of a phased array for 
the "red" channel relative to the "black" channel. It may also 
be achieved by use of a multiply-fed parabolic reflector, in 
which the unmodified beam patterns are used for the "black" 
channels, but in which three adjacent feeds are energized 
each with Vxd of the energy destined for a "red" cell. Due 
to coherent addition, the gain in the center of the "red* 1 cell 
will be 3 times a "black" beam gain at that point, effectively 
"filling in the hole". The "green" beams are formed in 
exactly the same way for the "green" frequency or timeslot. 
This is achieved using ground based hardware that directs 
the appropriate combination of signals through the transpon- 
der channel directly associated with each of the feeds. 

In FIG. 35, the smaller circles show the areas out to which 
a particular channel is to be used, beyond which a different 
channel is available with a more centrally directed beam. 
The area has been filled in in the case of the black beam to 
assist in identifying it. This area extends from the center of 
a beam out to the beam radius over root(3), that is the "cell" 
area is only Vsrd of the "beam" area, and mobiles in the 
"cell" only use the "beam" out to slightly more than 50% of 
the beam radius. 

In practice, of course, many more than three channels are 
available per beam, so it is possible to plan for cells that are 
only 1/M of the beam spot area where M is the number of 
channels available. If M=7 for example, beams are only used 
out to l/root(7) of their radius, as illustrated in FIG. 36. In 
practice M will be at least 100, so cell radius can be Vioth of 
the beam radius, hence the gain and C/I performance of the 
beam configuration is only important for a fraction of the 
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beam spot coverage. This does not necessarily mean it is 
possible to shrink the spot to obtain more gain, as it would 
not then be possible so easily to "fill in the holes". With a 
large number of offset beams it is desirable to phase the 

5 physical feeds to produce peak gain anywhere, and as 
indicated by FIGS. 24 and 30. The hardest place to obtain 
gain (by phasing only two feeds together) is midway 
between two spots, and the gain under those circumstances 
is maximized by choosing the beam edge crossover points as 

10 indicated on FIGS. 24 and 30. The gain between two beams 
is then twice the beam edge gain (e.g., 3 dB up) while the 
gain between three beams is three times the gain of one 
beam at that point. This explains the scaling used on FIGS. 
24 and 30 for comparing the gain at those three points. 

IS Thus, the error correction coding of rate between Yi and 
l A which is, in any case, contemplated for power efficiency 
reasons, can also provide tolerance of the C/Is obtained with 
immediate frequency re-use in every beam, if the technique 
of re-use partitioning just explained above is employed. The 

20 technique of re-use partitioning achieves this without resort 
to null-creation or interference cancellation, i.e., all antenna 
degrees of freedom are used to maximize gain. The tech- 
nique of interference cancellation or creating pattern nulls at 
the center of neighboring cells can be employed as a further 

25 bonus to reduce C/I from neighboring beams to negligible 
proportions. 

An exemplary coding scheme which can be used to 
implement this exemplary embodiment of the present inven- 
tion is punctured convolutional coding based on rate Yx or 

30 Vsth, but in which the coding of each uncoded speech bit is 
adapted according to its perceptual significance to a level 
between, for example, rate Yi and rate Vs. Although BPSK 
tolerates 3 dB lower C/I than QPSK, there seems no reason 
to incur the 2:1 bandwidth efficiency loss. The C/I tolerance 

35 of QPSK with twice the coding is in fact better than BPSK 
with half as much coding, therefore a quaternary modulation 
can be used at least for the downlink. 

The above discussion is based on coherent demodulation 
performance, which is achievable in the satellite-to-mobile 

40 channel using quite wideband TDM and is not achievable 
with narrowband FDM. The criterion for the downlink 
method is that the number of information bits to be demodu- 
lated and decoded shall be large over the time during which 
the fading component of the channel can be considered 

45 static, that is about 200 at 2.5 GHz and a vehicle speed 
of 100 Krn/Hr. Thus the information rate needs to be a 
couple of orders higher than 5 Kb/s and with, for example, 
a l A mean coding rate the transmitted bit rate needs to be 
greater than 1.5 Mb/s, which will pass through a 1 MHz 

50 bandwidth channel using quaternary modulation. The capac- 
ity provided by systems based on the foregoing techniques 
is of the order of 100 Erlangs per MHz per spot area using 
a 4 Kb/s vocoder or 166 Erlangs per MHz per spot using a 
2.4 Kb/s vocoder. 

55 One exemplary technique for implementing the above- 
described re-use scheme of the present invention is by 
means of ground-based beam-forming, as described earlier 
in this specification. This involves providing feeder links to 
carry the signal for each antenna feed from the central 

60 ground station to the satellite in such a way that the relative 
phase and amplitude differences between the signals is 
preserved. Using such a coherent transponder, only one 
transponder channel is needed on the satellite per antenna 
feed point. 

65 An alternative means of implementing the present inven- 
tion for the fixed-beam case disclosed is set forth below, 
which avoids the need for coherent feeder links at the 
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expense of more hardware on the satellite. The use of 
on-board beam forming is simplest if it is fixed and not 
variable, which may be more suited to a geostationary 
satellite that illuminates fixed areas. For a non-geostationary 
satellite this exemplary embodiment of the present invention 
can still be employed, but it is not then so easy to obtain the 
advantage of systematically adapting the beam forming to 
compensate for satellite motion so that beams illuminate 
fixed areas. 

The implementation of the fixed beam-forming transpon- 
der is shown in FIG. 37 for the FDMA case, i.e., the total 
available bandwidth is divided into N sub-bands, each of 
which is used to illuminate areas on the ground according to 
a cellular re-use pattern such as shown in FIG. 35. The case 
of three sub-bands — designated black, red and green as per 
FIG. 35 — is used by way of illustration. 

A set of transponder channels 37 receive signals from a 
corresponding set of feeder links and downconvert them to 
a suitable intermediate frequency for amplification and 
filtering. The I.F. outputs of 3710 are applied to I.F. beam- 
forming network 3720 which forms weighted combinations 
of the I.F. signals. The "black" channels are arbitrarily 
chosen to correspond directly to unmodified antenna 
patterns, i.e., black signal 1 shall be radiated directly and 
only through antenna feed number 1; black signal 2 shall be 
radiated only through antenna feed number 2, etc. The 
beam-forming network thus connects black channels with 
unity weighting into only that summing network corre- 
sponding to the designated antenna feed. 

The red channels and green channels however shall be 
radiated with a beam pattern centered midway between three 
black beams. The red beam that shall lie midway between 
black beams 1, 2 and 3 is thus connected to their associated 
three summing networks through voltage/current weightings 
of l/root(3) (power weightings of Vi). One third of the "red" 
energy is thus radiated through each of the three feeds 
surrounding the desired "red" center. Likewise, the green 
beam lying midway between black beams 2, 3 and 4 is 
connected via weightings of l/root(3) to the summers asso- 
ciated with feeds 2, 3 and 4. The weightings quoted above 
are exemplary and simplified for the purposes of illustration. 
Since the I.F. beamforming network can in principle be 
realized with a network comprising mainly simple resistive 
elements, more complex sets of weights can be used with 
acceptable complexity impact For example, a beam can be 
formed by feeding more than three adjacent feeds, and 
negative weights can be used to create nulls in the radiation 
pattern at desired places or otherwise to reduce the sidelobe 
levels in order to increase the C/I. 

One method to form a resistive I.F. beamforming network 
uses a continuous sheet or thin film of resistive material 
deposited on an insulating substrate. This sheet is notionally 
regarded as corresponding to the two dimensional surface to 
be illuminated by the beams. Signal currents corresponding 
to the "black" beam signals are injected into the sheet at 
points disposed in correspondence to the centers of the 
"black" cells, while "red" and "green" signal currents are 
injected at sets of points midway between the black signal 
injection points and each other, as per FIG. 35. FIG. 38 
illustrates the injection points by the labels T. 

Signal currents corresponding to the desired combinations 
of the black, red and green signals are extracted from the 
resistive plane by contacts disposed midway between the 
black, red and green injection points. These current extrac- 
tion points are indicated by 'E' in FIG. 38. This technique 
provides the same weight distributions for the black, red and 
green beams, in contrast with the previous example that had 
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single weights of 1 for the black beams and three equal 
weights of l/root(3) for the red and green beams. Extracted 
currents are fed to "virtual earth" amplifier inputs or low- 
input impedance amplifiers such as grounded base bipolar 

5 transistors. The set of weights realized by this technique can 
be tailored by choice of the shape and size of the current 
injection and extraction contact lands. No simple rule is 
proposed for deciding the size and shape — a proposition 
must simply be verified by carrying out a two-dimensional 

10 finite-element computer analysis of the current flow in and 
potentials existing on the resistive sheet. 

Once the combined signals have been produced by the I.F. 
beam forming network they are fed to a bank of upconver- 
tors 3730 for frequency translation to the desired satellite- 

15 mobile frequency band. The upconvertors are all driven by 
the same local oscillator signal so as to preserve the relative 
phase of the signals, and have matched gain to preserve 
relative amplitudes. The unconverted signals may then be 
amplified by a matrix power amplifier 3740 to raise the 

20 power level to the desired transmit power. 

The inventive technique described above can be extended 
to produce any number of virtual beams associated with 
subdivisions of the total frequency band available. In the 
three-color example, each "color" is associated with a '/3rd 

25 sub-bandwidth. If 16.5 MHz total is available, for example, 
each transponder channel bandwidth can be nominally 5.5 
MHz. If the number of feeds is 37, for example, 37, 5.5 MHz 
"black" beams are generated, 37, 5.5 MHz wide red beams 
and 37, 5.5 MHz green beams. Thus the total bandwidth 

30 available for communication is 37 times 16.5 MHz, as it 
would have been had it been possible to employ immediate 
frequency re-use of all 16.5 MHz in the "black" beams only. 
Thus the present invention provides the same efficient use of 
bandwidth as an immediate frequency re-use pattern but 

35 with a considerably improved C/I- 

The extra capacity thus available is, in the FDMA version, 
obtained by increasing the number of transponder channels 
and thus hardware complexity proportionally. It will now be 
shown how an exemplary TDMA embodiment can be 

40 advantageously constructed in which the capacity increase is 
obtained without increased hardware complexity. 

FIG. 39 illustrates the exemplary TDMA embodiment. In 
this case the number of transponder channels 3910 is the 
same as the number of antenna feeds, and the bandwidth of 

45 each channel is the full bandwidth available to the system. 
The I.F. beam forming network 3920 also functions as 
previously described to synthesize black, red and green 
beams, but only one color is connected at a time to the set 
of transponder channels by virtue of the commutating 

50 switches 3911. Either (1) all transponder channels are con- 
nected to a corresponding number of "black" beam inputs, 
or, (2) by operating switches 3911 all at the same time, all 
transponders are connected to the red beam inputs, or (3), as 
shown in FIG. 39, to the green beam inputs. 

55 The switches are cycled such that for a first portion of a 
TDMA frame period (e.g., VSrd) the black beams are used, 
for a second portion of the time the red beams are energized, 
and for a third portion of the time the green beams are 
energized. The time periods during which the switches dwell 

60 in each position do not have to be equal, and can in principle 
be adapted according to which color has the highest instan- 
taneous capacity demand in any cell. The functioning of the 
rest of the transponder is as previously described for the 
FDMA case. 

65 It will be appreciated that the commutation of switches 
3911 is synchronized with transmissions from the central 
ground station or stations, and this can be achieved by any 
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of a variety of techniques, such as providing an on-board 
clock that can be programmed form the ground to execute 
the regular cycle of switch operations and to synchronize the 
ground station transmissions to the satellite, which is the 
master timer. Alternatively a ground station can transmit a 
switch command using a control channel separate from the 
traffic channels. The method of achieving synchronism of 
the beam gyrations with the ground network is immaterial to 
the principle of the present invention. 

It may be appreciated that, although both the TDM A and 
FDMA versions of the invention disclosed above used fixed 
beam forming networks, it is possible by an obvious exten- 
sion of the method to permute the allocation of frequencies 
or time slots to beam colors by use of switches 3911 
controlled from the ground in such a way as to keep the areas 
of the earth illuminated by a given frequency or time slot as 
nearly as possible fixed. This is of course achieved more 
accurately by using greater numbers of "colors" (that is 
timeslots or sub-bands). Increasing the number of sub-bands 
involves hardware complexity in the FDMA case, so the 
TDMA version is preferred in this respect. The phasing of 
the commutation switches may thus be chosen so as to 
compensate for satellite motion and keep the areas illumi- 
nated by a particular timeslot or frequency more or less 
constant. The present invention can be applied to any 
number of timeslots and sub-bands, and in the latter case a 
digital implementation comprising analog-to-digital conver- 
sion of the transponder signals, digital filtering and digital 
beam forming using digital weight multiplication may be 
advantageous. 

The above-described exemplary embodiments are 
intended to be illustrative in all respects, rather than 
restrictive, of the present invention. Thus the present inven- 
tion is capable of many variations in detailed implementa- 
tion that can be derived from the description contained 
herein by a person skilled in the art. All such variations and 
modifications are considered to be within the scope and 
spirit of the present invention as defined by the following 
claims. 

What is claimed is: 
1. A method for communicating a plurality of first signals 
from a first station to a corresponding plurality of second 
stations and a plurality of second signals from said second 
stations to said first station involving the steps of: 

processing said first signals into a form representative of 45 

modulated radio signals; 
combining said representations to form a number of 
transmit signals using a set of combining parameters; 
transmitting said transmit signals using a transmitter and 

antenna for each of said signals; 
receiving said transmit signals at said second stations to 

produce received signals; 
measuring at said second stations an amount received of 
at least one unwanted signal relative to the amount of 55 
wanted signal received to produce measurement 
results; 

transmitting said measurement results from said second 
stations to said first station; 

receiving said measurement results at said first station and 
using said results to modify said combining parameters 
such that unwanted signals received at said second 
stations are reduced and wanted signals enhanced; 

wherein said step of measuring further includes the step of 
correlating said received signals with known signals to 
produce said amount of at least one unwanted signal, 
and said amount of wanted signal. 
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2. A method according to claim 1 in which said first 
station comprises a ground station communicating with at 
least one orbiting satellite. 

3. A method according to claim 2 in which said corre- 
sponding transmitters form part of a satellite transponder for 
relaying said transmit signals from said ground station via 
said satellite. 

4. A method according to claim 3 in which said transmit 
signals are relayed by a coherent transponder in such a way 
as to preserve their relative phases and amplitudes. 

5. A method according to claim 1, wherein said measuring 
step includes forming a ratio of said amount of at least one 
unwanted signal to said amount of wanted signal. 

6. A method according to claim 5, wherein said ratio 
forms a complex value. 

7. A method according to claim 1 in which said step of 
correlating said received signals with known signals com- 
prises correlating said received signals with a set of codes. 

8. A method according to claim 7, wherein said amount of 
wanted signal is produced by correlating said received 
signals with at least a first code, and said amount of at least 
one unwanted signal is produced by correlating said 
received signals with at least a second code. 

9. A method according to claim 1, wherein said step of. 
transmitting comprises transmitting only largest ones of said 
measurement results to said first station. 

10. A method according to claim 1, wherein said step of 
receiving at the first station modifies said combining param- 
eters such that wanted signals received at said second 
stations are maximized. 

11. A method for communicating a plurality of first signals 
from a first station to a corresponding plurality of second 
stations and a plurality of second signals from said second 
stations to said first station involving the steps of: 

processing said first signals into a form representative of 
modulated radio signals; 

combining said representations to form a number of 
transmit signals using a set of combining parameters; 

amplifying said transmit signals to desired transmit power 
levels and transmitting said transmit signals using a 
corresponding transmitter and antenna for each trans- 
mit signal; 

measuring at said second stations an amount received of 
at least one unwanted signal relative to the amount of 
wanted signal received to produce measurement 
results; 

transmitting said measurement results from said second 
stations to said first station; 

receiving said measurement results at said first station and 
using said results to modify said combining parameters 
such that unwanted signals received at said second 
stations are reduced and the sum of said transmit power 
levels needed to transmit said first signals is reduced. 

12. A method according to claim U in which said first 
station comprises a ground station communicating with at 
least one orbiting satellite. 

13. A method according to claim 12 in which said 
corresponding transmitters form part of a satellite transpon- 
der for relaying said transmit signals from said ground 
station via said satellite. 

14. A method according to claim 13 in which said transmit 
signals are relayed by a coherent transponder in such a way 
as to preserve their relative phases and amplitudes. 

15. A method of claim U, wherein said step of receiving 
said measurement results further reduces said sum of said 
transmit power levels to a minimum. 
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16. A system for wireless communication between at least wherein said logic for determining further includes logic 

one first station and at least one second station, comprising: for correlating said received signals with known signals 

a first station including a processor for processing first to produce said amount of at least one unwanted signal, 

signals into a form representative of modulated radio and said amount of wanted signal, 

signals; 5 17. A system according to claim 16, wherein said logic for 

a first station combiner for combining said representations determining includes forming a ratio of said amount of at 

to form a number of transmit signals using a set of least one unwanted signal to said amount of wanted signal, 

combining parameters; 18. A system according to claim 17, wherein said ratio 

a first station transmitter for transmitting said transmit forms a complex value. 

signals; 19. A system according to claim 16 in which said logic for 

a second 'station receiver for receiving said transmitted correlating said received signals with known signals corn- 
signals to produce received signals; logic for correlating said received signals with a set of 

second station logic for determining an amount received C °^n S ' A , t , . 1ft . j 

bysaidreceiverof atleastoneunwant ed signal relative is 20. A system according to claim 19, wherem said amount 

to the amount of wanted signal received by said of wanted signal is produced by correlating said received 

receiver to produce measurement results; s, S nals W1,h , a ' least a , first • nd . s V d am< T t ° f at lea ?! 

v . . one unwanted signal is produced by correlating said 

a second station transmitter for transmitting said mea- received with at least a second code. 

surement results from said second station to said first ^ 21 Asystem accordin g t0 claim 16, wherein said second 

station; station transmitter transmits only largest ones of said mea- 

a first station receiver for receiving said measurement surement results to said first station. 

results from said second station transmitter; 22. A system according to claim 16, wherein said first 

wherein said first station uses said results to modify said station modifies said combining parameters such that wanted 

combining parameters such that unwanted signals 25 signals received at said second station are maximized. 

received at said second stations are reduced and wanted 

signals enhanced; * * * * * 
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